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TEST MEASUREMENTS AND NASTRAN ANALYTICAL VALUES 
L. Purves, R. Strang, M. P. Dube, P. Alea, N. Ferragut, and D, Hershfeld 

ABSTRACT 

This document describes the software and procedures of a system of 
programs used to generate a report of the statistical correlation between 
NASTRAN modal analysis results and physical test results from modal surveys. 
Topics di r tussed include a brief mathematical description of statistical correla- 
tion, a user’s guide for generating a statistical correlation report, a programmer’s 
guide describing the organization and functions of individual programs leading 
to a statistical correlation report, and a set of examples including complete 
listings of programs, and input and output data. 
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STATISTICAL CORRELATION OE STRUCTURAL MODE SHAPES 


FROM 

TEST MEASUREMENTS AND NASTRAN ANALYTICAL VALUES 

Section 1 

INTRODUCTION AND OVERVIEW 

In recent years there has been a large improvement in the design of mechanical structures, 

The digital computer has significantly contributed to this improvement. Computer implementa- 
tions of the finite clement technique, in programs such as NASTRAN (NASA STRuctural ANalysis), 
have made it possible to predict accurately and in great detail the behavior of complex structures. 

Computer Aided Design (CAD) and Computer Aided Manufacturing (CAM) have lowered the 
costs of creating complex and sophisticated structures. Computer Aided Testing (CAT) has allowed 
much more extensive testing to be performed on structures. It has appeared that relatively less 
attention has been paid to the computerized comparisons between analytical and test data. Since 
many million bits of analytical and test data can be easily created with computer techniques, it is not 
practical to do complete comparisons manually. The use of appropriate computer methods, some of 
which can be based on statistical analyses, can increase the speed and completeness with which com- 
parisons between test and analytical data can be made. 

Potentially, automated techniques can help identify whether disagreement is due to faulty test 
procedures or faulty modeling. Another promising long range benefit of the automated correlation 
of analytical and test data is that structural design can be further improved. There is presently no 
widely accepted and effective way of validating analytical results on real problems. Statistical com- 
parisons with equivalent test data appear to provide an acceptable means for this validation. If this 
technique becomes widely used and improved, then it could gradually lead to improved analytical 
techniques, thus, leading in turn, to improved designs and eventually to improved structures. 
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The purpose of this document is to describe some recent developments, based on statistical 
correlation, of comparing two given sets of data, With the techniques described below, NASTRAN 
and several pre-processing programs place analytical results in a format in which they can be statis- 
tically compared with the test results, The STATCORR program has been developed to compile a 
statistical correlation report in which the engineer can easily interpret and make necessary decisions 
based on the statistical comparisons of the two sets of data. 

In aerospace structures, it is necessary to compare structural vibration properties obtained 
from testing with those obtained from analysis. The most meaningful comparisons of data arc those 
based on like vibrational modes, A test mode and an analytical mode are alike if the amplitude of 
the standing wave pattern is “exactly” the same at ail sample points, or more realistically, if the pair 
of modes correlate closely. 

The limited amount of experimental data and the limited access in which to install sensors at 
interior points are factors that constrain the points at which comparisons can be made. There is a 
disproportionately large number of points at which analytical data is available for comparison. The 
objective is to draw useful comparisons within the constraints of available data. The intent is to be 
able to identify the mode shapes from both of the sources; to analyze the relationship of any two 
mode shapes; to delineate at which frequencies the test mode and the analytical mode of a closely 
related pair appear; to give the locations at which two modes diverge beyond a threshold; and to 
isolate the greatest divergence. 

The chart in Figure 1 provides a pictorial overview of the NASTRAN software system and the 
statistical correlation software on the GSFC Code 750 VAX-1 1/780 computer. 

Section 2 of this document provides a brief mathematical description of the more essential 
quantities used in statistical analyses of data and in the STATCORR program. 
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Section 3 is u user’s guide covering the execution of NASTRAN and the statistical correlation 
utility programs. 

Section 4 is a programmer’s guide describing the organization and functions of individual pro- 
grams leading to the statistical correlation report. 

Section 5 contains a set of examples, including a description of the Spartan* 1 model. Also pro- 
vided are complete listings of some of the input and output data. 

The Appendix contains the NASTRAN input data of the SPARTAN-l model, and an example 
of a statistical correlation report. Also included are the FORTRAN source code listings of 
STATCORR and tire other preprocessing programs. 

The Statistical Correlation software is part of the NEXUS C'AD/CAM system at GSFC, and it is 
available from COSMIC. The files exists in the NliXUS ,LIB: [NEXUS.NASTRAN. STATCORR] 
directory, with the programs in the [.PROGRAMS] subdirectory and the demonstration files in the 
[.DEMO] subdirectory. The demonstrations are based on the SPARTAN- 1 model (described in 
detail in section 5), and can be run with the help of the STATCOR.R.COM command procedure file. 
The analytical NASTRAN input data files are called: FRE.NID, FXT.N1D, and ETUPLT.N1D; the 
experimental input data files arc called: FMS. . . 
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THEORY OF STATISTICAL CORRELATION 

Static and dynamic characteristics of mechanical structures arc determined both analytically 
and experimentally, In cither case, the results of such determinations arc subject to random varia* 
tion. For example, analytical modes of a structure will vary because of differences in interpretation 
of drawings, and representations of elastic relations and mass properties. Similarly, experimental 
results will vary because of instrumentation errors, data recording noise, and data processing tech- 
niques, Despite these variations, it is hoped that both the analytical and experimental results will 
tend to correlate. 

Analysis and testing results can be systematically compared, using well-established principles of 
statistics, particularly, the correlation coefficient. (Sec Freund, “Mathematical Statistics,") The 
criterion for deciding whet he; two mode shapes describe the same mode will be how close their 
correlation coefficient is to unity. 

Some of the more important quantities used in statistical analyses arc described below, The 
mean is defined as the first moment of a distribution about the origin. That is, 

* 4 t *< ' r(x i> = ir t *i <■> 

1° 1 1=1 

wb,rc for the purpose of this analysis, the quantity Xj is the amplitude of a vibration mode at a par- 
ticular point “i”. Function f(x) is the probability function for the amplitude xj. With the proba- 
bility that a vibration mode exists at each point at the time of the sample measurement, then the 
probability function reduces to one at each point, f(Xj) = 1. The symbol n is the total number of 
points in the sample where a measurement is made. 

The variance is defined as the second moment of a distribution about its mean. That is, 

* x = ~ £ (xj - x)'f(Xj) = - £ (x l-^ (2) 

1=1 1=1 
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For the purposes of this application, the RMS (root mean square) of u distribution Is defined as: 


RMS = standard deviation « SQRT(varlnncc) = o x (3) 


Covariance is defined as the first product moment of two separate distributions about tiicir re- 
spective means. That is, 

a xy = T 'll < x l “ *> * &i - y) W) 

JB 1 


A quantity closely related to the covariance is the correlation coefficient, p. It is defined as: 

- g xy 
p xy o x • o y 


(5) 


Tiie correlation coefficient is independent of any scaling of the mode shapes, and thus frees 
botli tiie analyst and the experimentalist to scale data in the manner of their ciioicc. Each set of 
data is scaled to Its own standard deviation and both arc placed on equal footing regardless of how 
tiie modal data was originally obtained. All modes in botli sets are compared with cacli other, and 
tiie correlation coefficient is computed for every combination. 

It is instructive to determine where and by how much mode shapes differ from one another. 
However, because of scaling, the difference between the amplitudes of two mode shapes is not sig- 
nificant. One method to obtain significant differences is to rescale the mode by their respective 
RMS values. This rescaling puts all mode shapes on an equal basis. The relative difference between 
the amplitudes of two mode shapes at a point is then 

Relative Difference = — • - — (6) 

p x Py 


Another method rescales one of two mode shapes to be compared so that the mean squared differ- 
ence between the two mode shapes is minimized. The scaling factor for this method is: 


( 7 ) 


C’ 


>' x,y| 

i° i 


x *> 

n « 1 


If the x mode shape Is multiplied by this sealing factor, the mean squared difference between the 
x and y mode shapes will equal: 

» a 




( 8 ) 


The square root of the mean squared difference is called the standard difterenec, S. The dillerence 
between the amplitudes of the two mode shapes can now be compared to the standard difference to 
detect any significant variations between mode shapes. The sealed difference is defined as: 


Scaled Difference = (C’Xj - yj)/S 


(9) 


Most of the above quantities are used in the statistical correlation report, and an example of 
such a report is given in Appendix C. The variances, correlation coefficients, and RMS values are 
computed in the CORRMS subroutine of the STATCORR program. The relative and scaled differ- 
ences are computed in the RMSDIFF subroutine. Further discussion on the theory of statistical cor- 
relation will be found in a future Technical Memoranduntfcallcd “Statistical Correlation Analysis 
For Comparing Vibration Data From Test and Analysis/’ by L R. Purves, D. J. Hershfetd, T. CJ. 
Butler, and R. F. Strong. 

*TM 86081 
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Section 3 


USER’S (JUII)E 


flic following is a guide to aid the user in the execution of NASTRAN and the statistical corre- 
lation utility programs. The STATOORR utility program generates a statistical correlation report 
from results of analytical NASTRAN data and experimental test data. Also needed as input to 
STATCORR is an external grid point list, that is, a set of grid points governing the location and orien- 
tation of the test sensors. 


The data that is entered as input into the STATC’ORR program must be entered according to a 
specific format. There are several preprocessing programs which are used to properly format this 
data, and they are described beiow. 


First, one needs to process the NASTRAN analytical data. The NASTRAN input data file 
usually contains an Executive Control Deck, a Case Control Deck, and also a Bulk Data Deck. This 
file may reside on a disk or a magnetic tape. The tape may contain data from a previous NASTRAN 
run, and if this is the case, then the necessary files are copied from the tape to the user’s disk. This 
tape should include the following files: a print file (EX I ,PRT;I , for example), a command proce- 
dure file (EX I .COM: 1 ), a “new problem tape” file (EX1 .PTP;I), and a NASTRAN input data file 
(EX1.NID;! ) The file, EX1.PTP;!, is the “new problem tape” of a previous NASTRAN run, but 
becomes an “old problem tape" for the next NASTRAN restart run. 


Figure 2 shows the command procedure for assigning files in a NASTRAN restart run. The 
default must first be set to the user’s directory in each of the following command procedure files. 
The first assignment is the NASTRAN print file, In the VAX system, data written to Logical Unit 6, 
FOR006, goes to the log Tile for the batch job, by default. The procedure assigns a filename to 
FOR006, and the NASTRAN print output is written to the assigned file rather than the batch job 
log file. In the next assignment, the NASTRAN old problem tape file is assigned to the Logical 
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Unit 8, FORQ08. In the last assignment, NASTRAN writes tl - Bulk Data from the previous 
NASTRAN run to the punch file, the Logical Unit 77, FOR077. The punch file becomes the Bulk 
Data Deck for the NASTRAN input data of the next NASTRAN run. (This is described in more 
detail below.) The last step in the command file results in the execution of the NASTRANX proce- 
dure. The “IiXliC" command is a global symbol defined for users of the Code 750 VAX-1 i/780 
computer system for executing commands in the DRAO:[USRLIB.COM] directory. The following 
definition is used: 

EXEC: = = “(" DRAO: [ USRLIB.COM ] COMC M 
The COMC.COM command procedure contains the following DCL command: 

S ("DRAO: [USRLIB.COM] Tl’ T2’ ’P3’ T4’ T5’ ’P6’ »P7’ T8’ 

NASTRANX is a command procedure which may contain up to five input parameters, PI , P2, 
P3, P4, and PS. If any particular parameter is to be used, all the preceding parameters must also be 
defined by the user, but following parameters may be ignored, PI is the user’s NASTRAN input file 
name; in this example it would be, EX1.NID. P3 is the maximum number of NASTRAN links to 
execute (default = 100,000,000). In order to generate the punch file only the first link is needed. 
Since P3 was specified, P2, the NASTRAN executable assignment file name 

(default = (" DRAO: [ USRLIB.COM] NSASSIGN), 
must also. By specifying P2 with the word “DUMMY”, the default is used. 

EX1.COM; 1 

$ SET DEF [ ] 

$ ASSIGN EX1.PRT FOR006 
S ASSIGN EX1.PTP FOR008 
S ASSIGN EX1.DIC FOR077 
$ EXEC NASTRANX EX1.NID DUMMY 1 
$ DEASSIGN FOR006 
S DEASSIGN FOR008 
$ DEASSIGN FOR077 
$ EXIT 


Figure 2. Example of a User’s Command Procedure for Assigning Files in a NASTRAN Restart Run 


l ; or user's of the GSI’C Code 750 VAX-1 1/780 computer system, the submit command for 
running batch jobs is given by: 

S SUBM1T/QUKUE=NASTRAN KX1.COM 

The NASTRAN input data file should contain a NASTRAN card, and it is the first card in the 
NASTRAN data deck. (An example of a NASTRAN input data deck is given in Appendix A.) The 
Restart deck follows the NASTRAN and ID cards, and then followed in order by the Executive 
Control Deck, the Case Control Deck, and the Bulk Data Deck. In this case, the Bulk Data is empty, 
but one will be generated from this run as a punch file. The only modification to be made to the 
NASTRAN input data file prior to the run is the ECHO card in the Case Control Deck, and this 
should be changed to: 

ECHO = PUNCH 

The job is then submitted to batch, and the resulting punch file, EX1.DIC, is appended to the 
NASTRAN input data, EX1.NID, as the NASTRAN Bulk Data Deck. The ENDDATA card is then 
shifted to the last card of the deck. The following modifications to NASTRAN Input Data file arc 
made: the NASTRAN card along with all of the Restart cards arc deleted; in the Executive Control 
Deck, the CHKPNT YES card is deleted, the DIAG 14 card is changed to DIAG 14,21,22, and the 
following alter cards are added: 

ALTER 41 
EXIT S 
ENDALTER 

In the Case Control Deck, the ECHO card is modified to read: 

ECHO - NONE 

The command procedure file is also modified, and Figure 3 shows how the command procedure file 
should read. Several print files are generated during the NASTRAN run, and the SEXEC CONCAT 
procedure concatenates different versions of files with the same filename, as in this case, 
EX1.PRT;!, EX1.PRT;2, 
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This job is then submitted to batch, and a concatenated print file is generated. A grid point 
list file must also be obtained or constructed, and it may be called GRDPT.LIS. This file is in a 
frec-field format, and it contains a list of all of the grid points in which test sensors were located 
during the test, Next, the preprocessing programs arc run, First, the TESETDMI program is run; it 
is designed to be run interactively, and it asks for the name of the NASTRAN print file and the name 
of the grid point list. A sample TESETDMI dialog is given as: 


S RUN TESETDMI 

FILENAME OF NASTRAN PRINT FILE WITH DIAG 21: EX1.PRT 
FILENAME OF MODAL SURVEY GRIDPOINT LIST: GRDPT.LIS 
FILENAME FOR TESET DMI BULK DATA RECORDS: TESET.DMI 
WARNING - NO INTERNAL DOF FOR GRID POINT 382-1 

WARNING - NO INTERNAL DOF FOR GRID POINT 382-2 

WARNING - NO INTERNAL DOF FOR GRID POINT 382-3 

WARNING - NO INTERNAL DOF FOR GRID POINT 407-1 

WARNING - NO INTERNAL DOF FOR GRID POINT 407-2 

WARNING - NO INTERNAL DOF FOR GRID POINT 407-3 

WARNING - NO INTERNAL DOF FOR GRID POINT 412-1 

WARNING - NO INTERNAL DOF FOR GRID POINT 412-2 

WARNING - NO INTERNAL DOF FOR GRID POINT 412-3 

FORTRAN STOP - PROCESSING COMPLETED 
S 


An output file, called TESET.DMI, is generated, and litis file is appended to the NASTRAN input 
data file. The file, TESET.DMI, contains the same list of grid points as the original file, GRDPT.LIS, 
but in the proper DMI format to be inserted into the NASTRAN Bulk Data Deck. The card 
ENDDATA is then shifted to the last card in the deck, and EX1.NID is modified: the DIAG 14,21,22 
card is changed to DIAG 14, and the alter cards changed as follows: 


ALTER 74 

TABPCM LAMA, , . ,//FM S 
ALTER 78 

PARTN P11IG, ,TESET/,PHITE, ./+1/1/2/2 S 
OUTPUT3 PHITE, TESET/ /-1/N l=VCT/N2=DOF S 
EXIT $ 

ENDALTER 


The command procedure file is also modified, and Figure 4 shows how the file should read. 

This job is then submitted to batch, and a concatenated print file and punch file are generated. How- 
ever, the user must first delete the old print files, otherwise they will be concatenated with the new 


ones. 
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EX1.C0M;1 


$ SET DEF [ ] 

$ ASSIGN EX1.PRT FOR006 
$ EXEC NASTRANX EX1.NID 
$ EXEC CONCAT EX1.PRT 
S DEASSIGN FOR006 
$ EXIT 


Figure 3. Example of a User’s Command Procedure for Assigning and Concatenating a Print File 


EX1.COM;3 

$ SET DEF [ ] 

$ ASSIGN EX1.PRT FOROQ6 
$ ASSIGN EX1.PCH FOR077 
$ EXEC NASTRANX EX1.NID 
$ EXEC CONCAT EX1.PRT 
$ EXEC CONCAT EX1.PCH 
S DEASSIGN FOR006 
$ DEASSIGN FOR077 
S EXIT 


Figure 4. Example of a User’s Command Procedure for 
Assigning and Concatenating a Print File and a Punch File 


The following programs are then run interactively. The UNPACKDMI program asks for the 
name of the concatenated punch file, HXl.PCH, and also the names of two matrices, TFSHT and 
PIUTE, A sample UNPACKDMI dialog is as follows: 

S RUN UNPACKDMI 

ENTER NAME OF NEXT DMI FILE: EX I, PCII 

ENTER MATRIX NAME. (‘“"ALL” FOR ALL MATRICES): TESET 

MATRIX FILE WRITTEN WITH MROWS, NCOLS: 

TESET ,MTX 1656 I 

DO YOU WANT ANOTHER MATRIX IN THIS FILE? (Y OR N): Y 
ENTER MATRIX NAME, ( U! "ALL’’ FOR ALL MATRICES): PIUTE 
MATRIX FILE WRITTEN WITH MROWS, NCOLS: 

PHITE ,MTX 135 56 

DO YOU WANT ANOTHER MATRIX IN THIS FILE? (Y OR N): N 

DO YOU WANT ANOTHER FILE? (Y OR N): N 

NO MORE FILES REQUESTED 

S S 

Two output files arc generated, TESET.MTX and PHITE. MTX. The file, PHITE.MTX, contains the 
mode shape displacements and it will later be an input file into the STATCORR program. The file, 
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TESET.MTX, is another file containing u list of test grid points, but in the wrong format. However, 
the GRDPTLST program is run and asks for the name of the grid point ID file, TESET.MTX, and 
generates the grid point list in the correct format. A sample GRDPTLST dialog is given as: 

S RUN GRDPTLST 

BNTKR GRID POINT ID MATRIX FILENAME: TESET.MTX 
ENTER GRID POINT ID OUTPUT LIST FILENAME: GRDPTF.LIS 
135 ENTRIES WRITTEN TO GRID POINT LIST 
FORTRAN STOP 
S 

The output file is called GRDPTF.LIS, and it will also be an input file to the STATCORR program. 

Finally the LAMA program is run, and it asks for the name of the punch file, and the name of the 

LAMA table. A sample LAMA dialog is given as: 

S RUN LAMA 

ENTER INPUT PUNCH FILE NAME: EX 1, PCM 
ENTER NAME OF LAMA TABLE: LAMA 
PLANE 1 SYMMETRY NOT SPECIFIED 
PLANE 2 SYMMETRY NOT SPECIFIED 
PLANE 3 SYMMETRY NOT SPECIFIED 

ENTER OUTPUT MATRIX FILENAME: LAMA.TBL 
FORTRAN STOP 
S 


An output file is generated, called LAMA.TBL, and this file contains the modal frequencies, masses, 
and stiffnesses, and it is another input entered into the STATCORR program. Therefore, the three 
files, LAMA.TBL, GRDPTF.LIS, and PHITE.MTX, contain all the necessary information (analytical 
NASTRAN data) in the proper format, and they are ail entered as input into the STATCORR 


program. 


Before the statistical correlation program can be run, the experimental data files are needed. 
These experimental data files usually exist on floppy disks, and the necessary files are copied to the 
user’s disk with the help of the FLX utility of the VAX-1 1 /7 SO computer system. 


Finally, the STATCORR program is run. This program is designed to run interactively, as with 
the other preprocessing utility programs, and a sample STATCORR dialog is given on page C-l of 
Appendix C. Ail the programs may also be run as a batch job with the help of a command proce- 
dure file. An example of a command procedure file used to run all the necessary programs from the 
initial input to the final output is shown in Figure 5. 
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STATCORR.COM; 1 

$ SET DEF [NASCAT] 

$ ASSIGN NEW.PRT FOR006 
S EXEC NASTRANX NEW.NID 
$ EXEC CONCAT NEW.PRT 
$ 

$ RUN (NASCAT.UTIUTY.STATCORR] TESETDMI 

NEW.PRT 

GRDPT.LIS 

TESET.DMI 

$ 

$ EDIT NEW.NID 
T "DIAG 14,21,22" 

DEL 

INSERT ; DIAG 14 
T "ALTER 41" 

DEL 

INSERT ; ALTER 74 

INSERT ; TABPCH LAMA, , , ,//FM $ 

INSERT ; ALTER 78 

INSERT ; PARTN Pl-ll G, , TESET/,PHITE, ./+1/1/2/2 $ 
INSERT ; 0UTPUT3 PHITE,TESET//-1/N1*VCT/N2*D0F $ 
T "ENDATA" 

DEL 

EXIT 

$ 

S EDIT ENDDATA.DAT 
INSERT ; ENDDATA 
EXIT 
$ 

$ APPEND ENDDATA.DAT TESET.DMI 
$ APPEND TESET.DMI NEW.NID 
$ 

$ DEL ENDDATA.DAT; 

$ 

$ ASSIGN NEW.PCH FOR077 
$ EXEC NASTRANX NEW.NID 
$ EXEC CONCAT NEW.PCH 
$ EXEC CONCAT NEW.PRT 

$ RUN [NASCAT.UTILITY.STATCORR] UNPACKDMI 

NEW.PCH 

TESET 

Y 

PHITE 

N 

N 

$ 


Figure 5. Example of a User’s Command Procedure to Run NASTRAN, the Pre-Processing 
Programs, and STATCORR as a Batch Job (Page 1 of 2) 
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$ RUN [NASCAT.UTILITY.STATCORRI LAMA 
NEW.PCH 
LAMA 
LAMA.TBL 
$ 

S RUN [NASCAT.UTILITY.STATCORRI GRDPTLST 

TESET.MTX 

GRDPTF.LIS 

S 

$ ASSIGN AC2.PRT SYSSOUTPUT 

$ RUN [NASCAT.UTILITY.STATCORRI STATCORR 

Y 

LAMA.TBL 

PHITE.MTX 

GRDPTF.LIS 

[NASCAT.EXPDATA] FMS001.AC2 
[NASCAT.EXPDATA] FMS002.AC2 
[NASCAT.EXPDATA] FMS003.AC2 
[NASCAT.EXPDATA] FMS004.AC2 
[NASCAT.EXPDATA] FMS005.AC2 
[NASCAT.EXPDATA] FMS006.AC2 
[NASCAT.EXPDATA] FMS007.AC2 
[NASCAT.EXPDATA] FMS008.AC2 
NONE 
N 
N 
N 

0.01 

S DELETE AC2.PRT; 

$ PRA AC2.PRT 
$ DEASSIGN SYSSOUTPUT 
S DEASSIGN FOR006 
$ DEASSIGN FOR077 
S EXIT 


Figure 5. (Page 2 of 2) 
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Section 4 


DESCRIPTION AND ORGANIZATION OF DATA 


The following section describes the organization and functions of individual data sets and 
programs leading to the statistical correlation report. Consideration has been given to analytical 
and test data coming from various sources. The STATCORR program is operated on data indepen- 
dent of its source, as long as it is entered according to a specific format. During this development, 
the analysis was done using NASTRAN, and the tests were performed by tiie CiSFC" Environmental 
Test and Integration Branch. 

At the outset, comparisons must be based on data gathered from identical sources. Liaison 
between test and analysis personnel prior to cither the generation of the NASTRAN model or the 
instrumentation of the test article must first establish the points and component directions at which 
data is to be sampled. This set of locations (referred to as the TESET vector) governs the location 
and orientation of test sensors as well as the entries in the analytical partitioning vector. 


Equally important as matching data points is matching boundary conditions, The most diffi- 
cult experimental boundary to enforce is clamped', the easiest to enact is free-free. Liaison between 
the test engineer and the analyst helps to uncover any disparity between the test set-up and the 
analytical boundary conditions. 


First, a list is obtained which contains the grid point ID’s and their associated DOF' components. 
From this list, an external gridpoint and DOF file is created by the user with all entries in list- 
directed ( free- field ) format. Entries arc inserted in grid. DOF pairs, one pair for each gridpoint 
and DOF-component which is instrumented in the modal survey. 


The internal sequencing of the grid points used by NASTRAN for matrix operations can bo 
different from the analyst’s original gridpoint sequencing. The partitioning vector, i.e., the THSHT 
vector, must be organized according to internal sequencing, so therefore a correspondence between 
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internal and external sequencing is necessary. This is accomplished by inserting DIAG 21 and 22 
cards in the Executive Control Deck, and an ALTLiR card to liXIT after the GP4 module, When 
DIAG 21 and 22 arc activated, the detailed correspondence between the internal and external 
sequencing for the various sets is printed. 

A program called TESETDMI.FOR will generate the TESET matrix in the proper DM I format, 
called TESET.DM1. Inputs to this program arc the external gridpoint list file in free-field format 
and the NASTRAN print file with the DIAG 2 1 output. TESET, DM1 is constructed in ordered 
pairs, the first member of which is a NASTRAN internal sequence number which corresponds to a 
component of a test point. The second member of the pair Is a Grid Point ID number belonging to 
the test point. TliSKT.DMI is then inserted into the NASTRAN Bulk Data for use in partitioning 
the PH1G matrix and obtain the PHITK matrix. The PillG matrix is the original total set of eigen- 
vectors of order G, which is to be partitioned (reduced), to the PIUTE matrix. The PIUTE matrix 
is the reduced set of eigenvectors pertaining only to the test set of grid points. 

In order to transfer the NASTRAN eigenvalue results into STATCORR, a set of DMAP 
ALTERS has to be included in the NASTRAN input deck before running the program. The DMAP 
module, TABPCI 1, delivers the modal frequencies, masses, and stiffnesses as part of the LAMA table 
in DTI format. The module PARTN docs the partitioning of the PHIG matrix by the TESET vector, 
and OUTPUT3 delivers the resulting partitioned mode shapes ■ s a single matrix called the PH1TE 
matrix. Both the LAMA table and the PllITE matrix together with the TESET vector are printed 
out in a punch file. 

Before inserting the LAMA table, the PIIITE matrix, and the external gridpoint list file into 
STATCORR, they must first be extracted from the punch file and put into the proper format. A 
set of three preprocessing programs docs this task. The LAMA program, with the name LAMA and 
the name of the punch file as input, searches the punch file for the LAMA tabl? in DTI format. 

Once found, the LAMA table is printed in a separate file, usually called LAMA.TBL. The 


UNPACKDMI program, with tlie names PHJTli and TIiSiiT, and the name of the punch file as input, 
searches the punch file for the matrices PiilTli ami TFSIiT in J)MI format, and prints them in two 
separate files, usually called PHITU.MTX and TESLT.MTX. The third program, CiRDPTLST.FOR. 
with TivSliT.MTX as input, prints in a proper format the external gridpofnt list in a file usually called 
GRDPTF.LIS. The three files, LAMA.TBL, PMITJi.MTX, and GRDPTF.LIS are the three input files 
containing all the necessary analytical NASTRAN data in proper formats to be used in STATCORR, 

The STATCORR program, as well as the other preprocessing programs, may be run inter* 
actively, or with the help of a user's command procedure file, they may also be run as a batch job, 
as shown in Figure 5. In order to compile and link all the necessary programs, one may use a com- 
mand procedure file called COMPLINK.COM; 1 , shown in Figure 6. 

COMPLINK.COM; 1 

$ FOR 'PV 

$ LINK 'PI' 

$ DEL 'PI', OBJ; 

S PURGE 'PI'.EXE 

Figure 6. Example of a User’s Command Procedure to Compile and Link a Computer Program 

There exists a number of options to be run with the STATCORR program. Tiie first is a 
request for a separate output print file. This print file, called the statistical correlation report, is 
printed by default to the system SYSSOUTPUT file. Therefore, in order to obtain a hard copy, the 
user must assign the SYSSOUTPUT file to a user’s disk file with the help of the ASSIGN command. 
An example of this is given in Figure 5 at the start of the STATCORR program. 

The first page of the statistical correlation report summarizes the interactive dialog. Part 2 on 
the following page contains a list of frequencies, masses, stiffnesses, and symmetries for the analyt- 
ical modes, which were obtained from the LAMA table. However, only the first 200 modes are 
processed and listed, along witii only tiie first 1000 DOFs. Next is a iist of frequencies, dampings, 
and symmetries for the experimental modes, Note that tire symmetry factors i e not considered in 
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this document, but they will be discussed in detail in a future T.M, (See bottom of page 2-3). With 
this set of data, STATC’ORR computes the correlation coefficients and RMS values for all possible 
analytical/experimental pairs. 

Part 3 of the statistical correlation report contains all the correlation coefficients for analytical 
versus experimental comparisons. Next, STATC’ORR determines the best experimental mode-shape 
match (highest correlation coefficient) for each analytical mode-shape, and also determines the best 
analytical mode-shape match for each experimental mode-shape. STATC’ORR also computes rela- 
tive deviations greater than threshold for each mode. The desired threshold Is a parameter which 
the user chooses with the default given by the value of 0.050. 

Part 4 contains the analytical mode shapes and their best experimental matches. Also listed are 
tlie relative deviations greater than threshold, and scaled differences greater than threshold. Part 5 
contains the experimental mode shapes and their l est analytical matches. 

The user may also request to have the analytical versus experimental symmetries to be con- 
sidered for each mode, P'inally, another option available to the user concerns the printing of the 
analytical and experimental mode-shape vectors. 

A flowchart showing the sequence of steps from the initial NASTRAN input data to the statis- 
tical correlation report is given in Figure 7. Source code listings for the five programs are given in 
Appendices D to 11. The NASTRAN input deck (with DIAG 21 and 22) for the SPARTAN-1 model 
(See Section 5 for a detailed description) is shown in Appendix A; the external grid point list file for 
this example is given in Figure 14; a partial listing of the NASTRAN print file (with DIAG 21 and 
22) is shown in Appendix B; a listing of the TESliT matrix in DM1 tormat is shown in Figure 15; a 
partial listing oTthc LAMA tabic is shown in Figure 16; the final external grid point listing, after it 
has been processed by the preprocessing programs, is shown in Figure 17; an example oT an input file 
containing experimental data is shown in Figure 18; and an example of a statistical correlation report 
is shown in Appendix C. 
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Section 5 

l-XAMPLliS THU SPARTAN-I MODFL 

In order to better understand the above sequence of steps in Figure 7, an example involving a 
realistic model will be given, The model, known as SPARTAN-1 , is described in detail below. The 
NASTRAN Input Data listing of SPARTAN- 1, with free boundary conditions and no mass mock- 
ups, is given in Appendix A, Included in the figures given below are complete listings of some of the 
input and output data to STATCORR and the other pre-processing programs. The Fortran source 
codes listings for these programs arc given in the Appendices, Also given is an example of a statisti- 
cal correlation report, Figure 8 shows a plot of the SPARTAN-1 model without the grid >oints. 



Figure 8, The Spartan-I Model 
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There exists a need to verify the accuracy of dynamic finite element models derived from 
structures scheduled to 11 y on the Space Shuttle, The SPARTAN-l payload is the first of the 
SPARTAN shuttle series. These SPARTAN payloads have evolved from the Sounding Rocket 
Experiments and can use the Shuttle environment to improve old experimental results or obtain 
new data not available with Sounding Rocket Technology. 

The payload consists of box-like structures in which half of the volume is occupied by support 
equipment: electronics, the attitude control system, the tape recorder, and batteries. The other half 
of the volume is occupied by the experiment, which is kinematically mounted to minimize distortion. 

A finite model of the structure was developed using NASTRAN. The model was built with 
enough details to simulate the mounting scheme of the experiment and batteries pack, and subse- 
quently allow stress anslysis of joint subassemblies. One of the requirements for payload flight cer- 
tification on the shuttle is a coupled-loads analysis using modal synthesis. It is important that the 
finite element model be verified to provide credibility to the load analysis results. 

The verification test extracted modal information from the structural test model in order to 
show the degree of correlation between the empirical and analytical models. This modal informa- 
tion includes natural frequencies, damping, and real mode shapes. 

The modal survey of the structural test model was performed in three configurations: 

Phase I, Flight SPARTAN frame with no mass mockups in a free boundary condition. 

Phase II. Flight SPARTAN frame with no mass mockups in a fixed boundary condition. 

Phase III. Engineering test unit with mass mockups in a fixed boundary condition. 

The above test sequence was chosen to minimize potential errors due to boundary condition 
effects. The free boundary condition specified in Phase I was achieved by suspending the test item 
from a facility crane with shock cord. The mass mockups used in Phase III were only to simulate 
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the approximate weights and locations of the flight subassemblies, The modal survey consisted of 
exciting the structure independently at three different points and measuring the response at thirty- 
live locations, The approximate location of the drive and response points is shown in Figure 9, 
Typical test setups for Phases 1, II, and 111 are shown in Figures 10, 1 1, and 12. respectively. Base- 
band frequency response functions (FRF) from (0-500 llz) were calculated from the measured ref- 
erence and response points. The modal information was then extracted. 


liaeh phase of the modal survey resulted in about 15-26 real inodes. A statistical correlation 
was performed on these mode shapes using the STATC’ORR program. As an example of the results, 
consider the first mode of vibration for each Phase (see Figure 13). It can be seen from the results 


shown in the figure that there is a high degree of correlation between the analytical and empirical 
mode shapes. Note that the correlation coefficient for Phases 11 and 111 is negative which indicates a 
180-degree phase reversal for these mode shapes. Also observe the fact that even though the varia- 
tions in the natural frequencies occur, such as a 20 r / frequency shift in the Phase 11 mode, the mode 


shapes are still well correlated. 


In general, for this series of tests, the lower-order modes correlate to a higher degree than the 
high-order modes. This could possibly be attributed to modal density observed and thus poor 
modal information. 


In the NHXUS,_L!B: [NF.XUS.NASTRAN.STATCORR.UFMOj directory of NUXUS, there 
exists the demonstrations files using the SPARTAN-1 model as the prime example. Included are 
NASTRAN input data files and experimental data files for the SPARTAN- 1 , along with the 
STATCORR.COM command procedure file to run the statistical correlation demonstration. 
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Figure 14. External Grid point Listing in Field-Free Format 
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Figure 15. Listing of the TESET Matrix in D.MI Format 


ORIGINAL PAGE H? 
OF POOR QUALITY 


LAMA . TBL ; 1 

LAMA 193 6 

O.OOOOOOOOE+OO 
O.OOOOOOOOE+OO 
O.OOOOOOOOE+OO 
O.OOOOOOOOE+OO 
O.OOOOOOOOE+OO 
O.OOOOOOOOE+OO 
4 . 096200 18E+01 
8 . 67981339E+01 
8 . 90685120E+01 
9 . 94205B56E+01 
1 . 27756042E+02 
1 . 34907944E+02 
1 . 42572708E+02 
1 . 50033951E+02 
1 . 51256454E+02 
1 . 62937897E+02 
1 . 66910095E+02 
1 . 71395691E+02 
1 . 7467053 2E+02 
2. 12637375E+02 
2. 15316406E+02 
2 . 25123245E+02 
2 . 30611511E+02 
2 . 32908356E+02 
2 , 33731552E+02 
2 . 39917892E+02 
2 . 44069534E+02 
2. 55830124E+02 
2 . 6595880 1E+0 2 
2 , 76783447E+02 
2 . 78470520E+02 
2 . 87401855E+02 
2 . 95211243E+02 
3. 08831482E+02 
3 . 22150085E+02 
3 . 32467621E+02 
3 . 40402924E+02 
3. 50125061E+02 
3 . 60949066E+02 
3 . 94879242E+02 
4 . 01737518E+02 


Figure 1 6. Partial Listing of the LAMA Table 
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original PAGE 13 
OF POOR QUALITY 


GRHPTP. LIS ; 1 


114-1 

1 

114-2 

2 

114-3 

3 

108-1 

19 

108-2 

20 

108-3 

21 

123-1 

97 

123-2 

98 

123-3 

99 

106-1 

109 

106-2 

110 

106-3 

111 

112-1 

115 

112-2 

116 

112-3 

117 

192-1 

127 

192-2 

128 

192-3 

129 

110-1 

139 

110-2 

140 

110-3 

141 

111-1 

145 

111-2 

146 

111-3 

147 

141-1 

205 

141-2 

206 

141-3 

207 

145-1 

211 

145-2 

212 

145-3 

213 

150-1 

367 

150-2 

368 

150-3 

369 

269-1 

385 

269-2 

386 

269-3 

387 

272-1 

391 

272-2 

392 

272-3 

393 

255-1 

397 

255-2 

398 

255-3 

399 

257-1 

409 

257-2 

410 


Figure 17. The Final External Gridpoint Listing to be 
Entered into the STATCORR Program (Page 1 of 3) 
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m r \ MAU Pft--*- 
if puOli quauhv 


264-1 

955 

264-2 

956 

264-3 

957 

4-1 

1027 

4-2 

1028 

4-3 

1029 

12-1 

1039 

12-2 

1040 

12-3 

1041 

9-1 

1069 

9-2 

1070 

9-3 

1071 

7-1 

1075 

7-2 

1076 

7-3 

1077 

97-1 

1189 

97-2 

1190 

97-3 

1191 

237-1 

1219 

237-2 

1220 

237-3 

1221 

94-1 

1231 

94-2 

1232 

94-3 

1233 

105-1 

1315 

105-2 

1316 

105-3 

1317 

40-1 

1333 

40-2 

1334 

40-3 

1335 

62-1 

1351 

62-2 

1352 

62-3 

1353 

102-1 

1447 

102-2 

1448 

102-3 

1449 

100-1 

1453 

100-2 

1454 

100-3 

1455 

44-1 

1549 

44-2 

1550 

44-3 

1551 

66-1 

1651 

66-2 

1652 

66-3 

1653 


Figure 17. (Page 2 of 3) 


a 


KCr- 

i 


i 

1 ORIGINAL PAGE tS 

OF POOR QUALITY 


257-3 

411 

181-1 

451 

181-2 

452 

181-3 

453 

177-1 

475 

177-2 

476 

177-3 

477 

148-1 

493 

148-2 

494 

148-3 

495 

154-1 

499 

154-2 

500 

154-3 

501 

153-1 

511 

153-2 

512 

153-3 

513 

152-1 

523 

152-2 

524 

152-3 

525 

231-1 

529 

231-2 

530 

231-3 

531 

244-1 

589 

244-2 

590 

244-3 

591 

234-1 

601 

234-2 

602 

234-3 

603 

155-1 

715 

155-2 

716 

155-3 

717 

1-1 

727 

1-2 

728 

1-3 

729 

245-1 

775 

245-2 

776 

245-3 

777 

212-1 

865 

212-2 

866 

212-3 

867 

216-1 

889 

216-2 

890 

216-3 

891 

261-1 

937 

261-2 

938 

261-3 

939 


Figure 17. (Page 3 of 3) 
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ORIGINAL P.VI* V:i 

of poor Q'ir.nv 


FMS001 . AC2 ? 1 


69.11099 
1. 1128000E-02 


0 

0 

236 

1 -3. 

234 

1 -3. 

245 

1 -6. 

244 

1 -5. 

257 

1 -2. 

255 

1 -2. 


0 

6037002E-02 

8017001E-02 

5154999E-02 

7073999E-02 

9363001E-02 

5717000E-02 


1 iguie 18. Example ol an Input File Containing Experimental D; 


ORIGINAL PAG? J3 

OF POOR QUALITY 




APPENDIX A. LISTING OF THE NASTRAN INPUT DATA (WITH DIAG 21 AND 22) 

FOR THE SPARTAN- 1 MODEL 


ID DSPPLT,FREE 
APP DISP 
SOL 3,0 
TIME 1500 
DIAG 14,21,22 
ALTER 41 
EXIT $ 

ENDALTER 

CEND 

TITLE DSP ANALYSIS — 

SUBTITLE = FREE-FREE MODES 

SPC = 1500 
METHOD = 100 
ECHO = NONE 
OUTPUT 

MAXLINES « 1000000 
LINE = 60 
VECTOR = ALL 


BEGIN 

BULK 









ASET1 

1 

234 

236 

244 

245 

255 

257 



ASET1 

2 

156 

158 

174 

175 

176 

183 

185 


ASET1 

2 

177 

179 







ASET1 

2 

199 

201 

217 

218 

219 

226 

228 


ASET1 

2 

214 

216 







ASET1 

3 

40 

42 

44 






ASET1 

3 

62 

64 

66 






ASET1 

3 

119 

121 







ASET1 

3 

137 

139 







ASET1 

123 

1 

4 

7 

9 

12 

15 

19 


ASET1 

123 

21 

28 

37 

47 





ASET1 

123 

51 

55 

59 

69 





ASET1 

123 

73 

84 

88 

90 

94 

97 

100 


ASET1 

123 

102 

105 

106 

108 

110 

111 



ASET1 

123 

112 

116 

117 

123 

125 




ASET1 

123 

132 

135 

141 

143 

147 




ASET1 

123 

148 

150 

152 

153 

154 

173 



ASET1 

123 

181 

220 

261 

264 





ASET1 

123 

269 

272 

376 

377 





CBAR 

1 

101 

1 

2 

0.0 

- 1 . 

1. 

1 

+1 

+1 



0 . 

-.6212 

.6212 

0 . 

-.6212 

.6212 


CBAR 

2 

101 

2 

3 

0.0 

- 1 . 

1 . 

1 

+2 

+2 



0 . 

-.6212 

.6212 

0 . 

-.6212 

.6212 


CBAR 

3 

101 

3 

4 

0.0 

— 1 . 

1 . 

1 

+3 

+3 



0 . 

-.6212 

.6212 

0 . 

-.6212 

.6212 


CBAR 

4 

101 

4 

5 

0.0 

- 1 . 

1 . 

1 

+4 

+4 



0 . 

-.6212 

.6212 

0 . 

-.6212 

.6212 


CBAR 

5 

101 

5 

6 

0.0 

- 1 . 

1 . 

1 

+5 

+5 



0 . 

-.6212 

.6212 

0 . 

-.6212 

.6212 



A~1 




(a 


CBAR 

+6 

CBAR 

+7 

CBAR 

+8 

CBAR 

+9 

CBAR 
+10 
CBAR 
+11 
CBAR 
+12 
CBAR 
+13 
CBAR 
+14 
CBAR 
+15 
CBAR 
+16 
CBAR 
+17 
CBAR 
+18 
CBAR 
+19 
CBAR 
+20 
CBAR 
+21 
CBAR 
+22 
CBAR 
+23 
CBAR 
+24 
CBAR 
+25 
CBAR 
+26 
CBAR 
+27 
CBAR 
+28 
CBAR 
+29 
CBAR 
+30 
CBAR 
+31 


29 

30 

31 


original i-vu-i? r? 

OF POOR QUALITY 


6 

101 

6 

7 

101 

0 . 

7 

8 

101 

0 . 

8 

9 

101 

0 . 

9 

10 

101 

0 , 

10 

11 

101 

0 . 

11 

12 

101 

0 . 

21 

13 

101 

.6212 

36 

14 

101 

.6212 

47 

15 

101 

.6212 

58 

16 

101 

.6212 

69 

17 

101 

.6212 

81 

18 

101 

• 6212 
105 

19 

101 

0 . 

104 

20 

101 

0 . 

103 

21 

101 

0 . 

102 

22 

101 

0 . 

101 

23 

101 

0 . 

100 

24 

101 

0 . 

99 

25 

101 

0 . 

98 

26 

101 

0. 

97 

27 

101 

0 . 

96 

28 

101 

0 . 

95 


7 

8 
9 



0.0 

-1 , 

6212 

.6212 

0 . 


0.0 

- 1 . 

6212 

.6212 

0 . ' 


0.0 

- 1 . 

6212 

.6212 

+ t 
0 . 


10 
-.6212 
11 


0.0 
.6212 


0 . 


12 


101 

101 

101 


0 . 

84 

-.6212 

70 

-.6212 

59 

'•6212 


36 

0 . 

47 

0 . 

58 

0 . 

69 

0 . 

81 

0 . 

90 

0 . 

104 
.6212 
103 
.6212 
102 
.6212 
101 
• 6212 
100 
.6212 
99 

.6212 

98 

.6212 

97 

.6212 

96 

.6212 

95 

.6212 

94 

.6212 

70 


6212 

u.o 

.6212 

- 1 . 

0 . 


0.0 

-1 

6212 

• 6212 

4 . # 
0 . 


1 . 

0.0 


.6212 

.6212 


1 . 

0,0 


.6212 

.6212 


1 . 

-.6212 

1 . 

-.6212 

1 . 

-.6212 

1 . 

-.6212 

1 . 

-.6212 

1 . 

-.6212 

1 . 


0 . 

59 

0 . 

48 

0 . 


1 . 
.6212 
1 . 

.6212 

1 . 

.6212 

1 . 

.6212 
0.0 
.6212 
0.0 
.6212 
0.0 
.6212 
0.0 
.6212 
0.0 
.6212 
0.0 
.6212 
0.0 
• 6212 
0.0 
.6212 
0.0 
.6212 
0.0 
.6212 
0.0 
• 6212 
- 1 . 

.6212 

- 1 . 

.6212 

- 1 . 

.6212 


0.0 
• 6212 
0.0 
.6212 
0.0 
.6212 
0.0 
.6212 
1 . 


0 . 

1 . 

0 . 

1 . 

0 . 

1 . 

0 . 

1 . 

0 . 

1 . 

0 . 

1 . 

0 . 

1 . 

0 . 

1 . 

0 . 

1 . 

0 . 

1 . 

0 . 

0.0 

-.6212 

0.0 

-.6212 

0.0 

-.6212 


0 . 

1 . 

0 . 

1 . 

0 . 

1 . 

0 . 

1 . 

0 . 

1 . 

0 . 

1 . 

.6212 

1 . 

.6212 

1 . 

.6212 

1 . 

.6212 

1 . 

.6212 

1 . 

.6212 

1 . 

.6212 

1 . 

.6212 

1 . 

.6212 

1 . 

■ 6212 
1 . 

.6212 

1 . 


0 . 

1 . 

0 . 

1 . 

0 . 


1 

.6212 
1 

.6212 
1 

.6212 
1 

• 6212 
1 

.6212 
1 

.6212 
1 

.6212 

1 

.6212 

1 

.6212 

1 

.6212 

1 

.6212 

1 

.6212 

1 

.6212 

1 

.6212 

1 

.6212 

1 

.6212 

1 

.6212 
1 

.6212 
1 

.6212 
1 

.6212 
1 

.6212 
1 

.6212 
1 

.6212 
1 

.6212 
1 

.6212 
1 

.6212 


+6 
+ 7 
+8 
+9 
+10 
+11 
+12 
+13 
+14 
+15 
+16 
+17 
+18 
+19 
+20 
+21 
+22 
+23 
+24 
+25 
+26 
+27 
+28 
+29 
+30 
+31 


A-2 




B iPTj 




ORIGINAL P/iGS 13 ’ 
OF POOR QUALITY 


CBAR 

+32 

CBAR 

+33 

CBAR 

+34 

CBAR 

+35 

CBAR 

+36 

CBAR 

+37 

CBAR 

+38 

CBAR 

+39 

CBAR 

+40 

CBAR 

+41 

CBAR 

+42 

CBAR 

+43 

CBAR 

+44 

CBAR 

+45 

CBAR 

+46 

CBAR 

+47 

CBAR 

+48 

CBAR 

+49 

CBAR 

+50 

CBAR 

+51 

CBAR 

+52 

CBAR 

+53 

CBAR 

+54 

CBAR 

+55 

CBAR 

+56 

CBAR 

+57 


32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 


101 

101 

101 

102 

102 

102 

102 

102 

102 

102 

102 

102 

102 

102 

102 

102 

102 

102 

102 

102 

102 

102 

102 

102 

102 

102 


48 37 

-.6212 0. 
37 25 

-.6212 0. 
15 2 

-.6212 0. 

37 38 

0. 0, 

38 39 

0. 0. 

39 40 

0 . 0 . 

40 41 

0, 0, 

41 42 

0* 0. 

42 43 

0 . 0 . 

43 44 

0. 0. 

44 45 

0 . 0 . 

45 46 

0 . 0 . 

46 47 

0 . 0 . 

59 60 

0. 0. 

60 61 

0. 0. 

61 62 

0- 0. 

62 63 

0. 0. 

63 64 

0 . 0 . 

64 65 

0. 0. 

65 66 

0. 0. 

66 67 

0 . 0 . 

67 68 

0 . 0 . 

68 69 

0 . 0 . 

4 16 

o. 0. 

16 28 

0. 0. 

28 40 

0- 0. 


-1, 
.6212 
- 1 . 
.6212 
-1. 
.6212 
0.0 
.5906 
0.0 
.5906 
0.0 
.5906 
0.0 
.5906 
0.0 
.5906 
0.0 
.5906 
0.0 
.5906 
0.0 
.5906 
0.0 
.5906 
0.0 
,5906 
0.0 
.5906 
0.0 
.5906 
0.0 
.5906 
0.0 
.5906 
0.0 
.5906 
0.0 
.5906 
0.0 
.5906 
0.0 
.5906 
0.0 
,5906 
0.0 
.5906 
0.0 
.5906 
0.0 
.5906 
0.0 
.5906 


0.0 

-.6212 

0.0 

-.6212 

0.0 

-.6212 

0.0 


0 . 

0.0 

0 . 

0.0 

0 . 

0.0 

0 . 

0.0 

0 , 

0.0 

0. 

0.0 

0 . 

0.0 

0 . 

0.0 

0 . 

0.0 

0 . 

0.0 

0 . 

0.0 

0 . 

0.0 

0 . 

0,0 

0 . 

0.0 

0 . 

0.0 

0 . 

0.0 

0 . 

0.0 

0 . 

0.0 

0 . 

0.0 

0 . 


0.0 

0. 

0,0 

0. 

0.0 

0. 


1. 
0. 
1. 
0 . 
1. 
0 . 
1 . 
0 . 
1. 
0 . 
1 . 
0. 
1 . 
0 . 
1. 
0 . 
1 . 
0 . 
1 . 
0. 
1 . 
0 . 
1 . 

0 , 

1, 

0. 

1 . 

0. 

1 . 

0. 

1 . 

0. 

1 . 

0 . 

1 . 

0. 

1. 

0 . 

1 . 

0. 

1. 

0. 

1 . 

0. 

1 . 

0 . 

1. 

0. 

1 . 

0 . 

1 . 

0 . 


1 +32 

.6212 

1 +33 

.6212 

1 +34 

■ 6212 

1 +35 

.5906 

1 +36 

.5906 

1 +37 

.5906 

1 +38 

.5906 

1 +39 

.5906 

1 +40 

.5906 

1 +41 

.5906 

1 +42 

.5906 

1 +43 

. 5906 

1 +44 

.5906 

1 +45 

.5906 

1 +46 

.5906 

1 +47 

.5906 

1 +48 

.5906 

1 +49 

.5906 

1 +50 

.5906 

1 +51 

.5906 

1 +52 

.5906 

1 +53 

.5906 

1 +54 

.5906 

1 +55 

.5906 

1 +56 

.5906 

1 +57 

.5906 


A-3 


ORIGINAL PAGE 
OF POOR QUAUTr 


£ 

t- 

CBAR 5Q 

+58 

102 

40 

a* 

; * 
fcr 

CBAR 59 

+59 

102 

0. 

51 

a 

CBAR 60 

+60 

102 

0. 

62 

< * 

CBAR 61 

+61 

102 

0. 

73 

% 
■ r 

CBAR 62 

+62 

102 

0. 

85 

V 

CBAR 63 

+63 

102 

0. 

7 

(t 

CBAR 64 

+64 

102 

0. 

19 

t 

CBAR 65 

+65 

102 

0. 

31 

**■ 

CBAR 66 

+66 

102 

0. 

42 

1 

■» 

CBAR 67 

+67 

102 

0. 

53 


CBAR 60 

+68 

102 

0. 

64 

ii* 

5 

CBAR eg 

+69 

102 

0, 

76 

* 

r 

CBAR 70 

+70 

102 

0. 

88 

\ 

CBAR 7i 

+71 

102 

0. 

9 


CBAR 72 

+72 

102 

0. 

22 


CBAR 73 

+73 

102 

0. 

33 


CBAR 74 

+74 

102 

0. 

44 


CBAR 75 

+75 

102 

0, 

55 


CBAR 76 

+76 

102 

0. 

66 


CBAR 77 

+77 

102 

0. 

78 


CBAR 7Q 

+78 

102 

0. 

91 


CBAR 79 

+79 

101 

0. 

25 


CBAR go 

+80 

103 

-.6212 

15 


CBAR 02 

+82 

103 

0. 

16 


CBAR 84 

+84 

103 

0. 

19 


CBAR q 5 

+85 

103 

0. 

22 


0. 




51 

0, 

62 

0 . 

73 

0. 

85 

0. 

97 

0. 

19 

0. 

31 

0. 


42 

0 . 

53 

0 . 

64 

0 . 

76 

0. 

8a 

o. 

100 

o. 

22 

0. 

33 

0. 

44 

0 . 

55 

0 . 

66 

0. 

78 

0 . 

91 

0 . 

102 

0 . 

15 
0 , 

16 
0 . 

19 

0 . 

22 

0 . 

21 

0 . 


* 

0.0 


.5906 
0.0 
.5906 
0.0 
.5906 
0.0 
.5906 
0.0 
.5906 
0.0 
.5906 
0.0 
.5906 
0.0 
.5906 
0.0 
.5906 
0.0 
.5906 
0,0 
.5906 
0.0 
.5906 
0.0 
.5906 
0.0 
.5906 
0.0 
.5906 
0.0 
.5906 
0.0 
.5906 
0.0 
.5906 
0.0 
.5906 
0.0 
.5906 
- 1.0 
.6212 
0.0 

-1.627 
0.0 

-1.627 
0. 0 

-1.627 

0.0 

-1.627 


O.l 
0 . 
0 .( 
0 . 
O.C 
0 , 
O.o 
0 . 
O.o 
0 . 
o.o 
0 . 
0,0 
0 . 

0 , o 
0. 
O.o 
0 . 
o.o 
0 . 
0.0 
0 . 
o.o 
0 . 

0.0 

0 , 

0.0 

0 . 

0.0 

0 . 

o.o 

0 . 

0.0 

0 . 

o.o 

0 . 

o.o 

0 . 

o.o 

0 . 

o.o 


0.1 

0. 

o.< 



1 +58 

.5906 

1 +e:q 


1 

.5906 

+59 

1 

.5906 

+60 

1 

• 5906 

+61 

1 

.5906 

+62 

1 

.5906 

+63 

1 

.5906 

+64 

1 

.5906 

+65 

1 

.5906 

+66 

1 

.5906 

+67 

1 

.5906 

+68 

1 

.5906 

+69 

1 

.5906 

+70 

1 

.5906 

+71 


1 

.5906 
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68 

257 

69 

258 

70 

SEQGP 

261 

157 

265 

125 

267 

93 

269 

65 

SEQGP 

263 

128 

270 

37 

271 

39 

376 

85 

SEQGP 

268 

94 

272 

66 

274 

40 

262 

124 

SEQGP 

273 

38 

260 

195 

264 

160 

266 

129 

SEQGP 

377 

113 

26 

156 

49 

159 

71 

168 

SEQGP 

1191 

5 

1232 

63 

251 

95 

254 

43 

SPCI 

1500 

4 

234 

235 

236 




SPC1 

1500 

4 

255 

256 

257 




SUPORT 

12 

123456 








ENDDATA 




APPENDIX B. Partial NASTRAN Print File Liotlng with DIAG 21 and 22 


— DSP ANALYSIS - AUGUST 24, 1983 

FREE-FREE MODES 

AAA USER INFORMATION MESSAGE 2118, SUBROUTINE GP4PRT - 
DIAG 21 SET -DOF VS. DISP SETS FOLLOWS. 

INT DOF EXT GP. DOF SBSG L A F N G ROS 


• • * • 


1620 

77 

- 

6 

101 



1620 

1620 

1 

1621 

54 

- 

1 



1514 

1621 

1621 

1325 

1622 

54 

- 

2 



1515 

1622 

1622 

1326 

1623 

54 

- 

3 



1516 

1623 

1623 

1327 

1624 

54 

- 

4 





1624 


1625 

54 

- 

5 





1625 


1626 

54 

- 

6 

102 



1624 

1626 

1 

1627 

55 

- 

1 

104 

190 

1517 

1625 

1627 


1628 

55 

- 

2 

105 

191 

1518 

1626 

1628 


1629 

55 

- 

3 

106 

192 

1519 

1627 

1629 


1630 

55 

- 

4 



1520 

1620 

1630 

1328 

1631 

55 

- 

5 



1521 

1629 

1631 

1329 

1632 

55 

- 

6 



1522 

1630 

1632 

1330 

1633 

70 

- 

1 



1523 

1631 

1633 

1331 

1634 

78 

- 

2 



1524 

1632 

1634 

1332 

1635 

78 

- 

3 



1525 

1633 

1635 

1333 

1636 

78 

- 

4 



1526 

1634 

1636 

1334 

1637 

70 

- 

5 



1527 

1635 

1637 

1335 

1630 

78 

- 

6 



1528 

1636 

1638 

1336 

1639 

67 

- 

1 



1529 

1637 

1639 

1337 

1640 

67 

- 

2 



1530 

1638 

1640 

1338 

1641 

67 

- 

3 



1531 

1639 

1641 

1339 

1642 

67 

- 

4 



1532 

1640 

1642 

1340 

1643 

67 

- 

5 



1533 

1641 

1643 

1341 

1644 

67 

- 

6 



1534 

1642 

1644 

1342 

1645 

65 

~ 

1 



1535 

1643 

1645 

1343 

1646 

65 

- 

2 



1536 

1644 

1646 

1344 

1647 

65 

- 

3 



1537 

1645 

1647 

1345 

1648 

65 

- 

4 



1538 

1646 

1648 

1346 

1649 

65 

- 

5 



1539 

1647 

1649 

1347 

1650 

65 

- 

6 



1540 

1640 

1650 

1348 

1651 

66 

- 

1 



1541 

1649 

1651 

1349 

1652 

66 

- 

2 



1542 

1650 

1652 

1350 

1653 

66 

- 

3 

187 

193 

1543 

1651 

1653 


1654 

66 

- 

4 



1544 

1652 

1654 

1351 

1655 

66 

- 

5 



1545 

1653 

1655 

1352 

1656 

66 

- 

6 



1546 

1654 

1656 

1353 

— COLUMN 

TOTALS 

--- 6 

102 187 

193 

1546 

1654 

1656 

6 1353 


V** 

bsf 

w 

3 

0 | 

& 

i 


B-2 





u 




/v" 


: -i- f 


- DSP ANALYSIS MGU3T 2,1 ' 1903 

free-free modes 


- 1 - - 2 - 
1„ 34-4 S4-5 



-1- 

-2- 

1° 

113-6 

1191-5 

11“ 

189-5 

183-5 

21 = 

257-4 

258-4 

31= 

238-4 

251-4 

41 = 

169-5 

167-5 

51 = 

246-4 

219-5 

61= 

217-5 

13-6 

71 = 

71-6 

200-5 

Bl = 

82-6 

237-4 

91 = 

86-6 

74-6 

101 = 

92-6 

89-6 


MFC DISPLACEMENT SET 
-3- -4- -5- 


SPC 

DISPLACEMENT SET 

-3- 

-4- 

-5- 

115-6 

193-5 

191-5 

187-5 

185-5 

184-5 

249-4 

250-4 

243-4 

252-4 

253-4 

244-4 

228-5 

222-5 

226-5 

240-4 

241-4 

235-4 

14-6 

26-6 

49-6 

204-5 

200-5 

206-5 

17-6 

18-6 

24-6 

87-6 

30-6 

35-6 

57-6 

80-6 

79-6 


-6- -7- -8- 


-6- 

-7- 

-8- 

127-6 

133-6 

129-6 

144-6 

146-6 

1232-5 

174-5 

175-5 

176-5 

239-4 

234-4 

165-5 

227-5 

224-5 

156-5 

157-5 

150-5 

162-5 

209-5 

210-5 

201-5 

27-6 

50-6 

83-6 

23-6 

20-6 

195-5 

34-6 

32-6 

93-6 

56-6 

77-6 

54-6 


B-3 


UY 


- DSP AHAVJTiS — AUGUST 24, 1983 

FREE-FREE 40DE3 


OM^t’ DISPLACEMENT SET 



-1- 

-2- 

-3- 

-4- 

-5- 

-6- 

-7- 

-8- 

1= 

114-1 

114-2 

114-3 

114-4 

114-5 

114-6 

113-1 

113-2 

11= 

113-5 

107-1 

107-2 

107-3 

107-4 

107-5 

107-6 

108-4 * 

21 = 

1191-1 

1191-2 

1391-3 

1191-4 

1191-6 

115-1 

115-2 

115-3 

31 = 

109-1 

109-2 

109-3 

109-4 

109-5 

109-6 

116-4 

116-5 

41 = 

118-2 

118-3 

118-4 

118-5 

118-6 

117-4 

117-5 

117-6 

51 = 

119-4 

119-5 

119-6 

128-1 

128-2 

128-3 

128-4 

128-5 

61 = 

120-2 

120-3 

120-4 

120-5 

120-6 

121-1 

121-2 

121-4 

71 = 

122-1 

122-2 

122-3 

122-4 

122-5 

122-6 

132-4 

132-5 

01 = 

123-5 

123-6 

124-1 

124-2 

124-3 

124-4 

124-5 

124-6 

91 = 

106-6 

112-4 

112-5 

112-6 

193-1 

193-2 

193-3 

193-4 

101= 

192-2 

192-3 

192-4 

192-5 

192-6 

191-1 

191-2 

191-3 

111 = 

110-4 

110-5 

110-6 

111-4 

111-5 

111-6 

127-1 

127-2 

121 = 

127-5 

136-1 

136-2 

136-3 

136-4 

136-5 

136-6 

126-1 

131 = 

126-4 

126-5 

126-6 

133-1 

133-2 

133-3 

133-4 

133-5 

141 = 

137-4 

137-5 

137-6 

129-1 

129-2 

129-3 

129-4 

129-5 

151 = 

138-3 

138-4 

130-5 

138-6 

130-1 

130-2 

130-3 

130-4 

161 = 

131-1 

131-2 

131-3 

131-4 

131-5 

141=4 

141-5 

141-6 

171= 

145-3 

145-4 

145-5 

145-6 

270-1 

270-2 

270-3 

270-4 

181 = 

273-1 

273-2 

273-3 

273-4 

273-5 

273-6 

271-1 

271-2 

191 = 

271-5 

271-6 

274-1 

274-2 

274-3 

274-4 

274-5 

274-6 

201= 

125-6 

134-1 

134-2 

134-3 

134-4 

134-5 

134-6 

254-1 

211= 

204-5 

254-6 

182-1 

102-2 

182-3 

182-4 

182-5 

182-6 

221 = 

190-3 

190-4 

190-5 

190-6 

189-1 

189-2 

189-3 

109-4 

231 = 

183-3 

183-4 

183-6 

108-1 

188-2 

188-3 

188-4 

188-5 

241= 

187-2 

187-3 

107-4 

187-6 

186-1 

186-2 

186-3 

106-4 

251 = 

185-1 

185-3 

185-4 

185-6 

184-1 

184-2 

184-3 

184-4 

261 = 

144-2 

144-3 

144-4 

144-5 

135-4 

135-5 

135-6 

143-4 

271 = 

M 2-1 

142-2 

142-3 

142-4 

142-5 

142-6 

149-1 

149-2 

281 = 

149-5 

149-6 

146-1 

146-2 

146-3 

146-4 

146-5 

139-1 

291 = 

139-5 

139-6 

147-4 

147-5 

147-6 

140-1 

140-2 

140-3 

301 = 

140-6 

150-4 

150-5 

150-6 

1232-1 

1232-2 

1232-3 

1232-4 

311 = 

151-2 

151-3 

151-4 

151-5 

151-6 

269-4 

269-5 

269-6 

321 = 

272-6 

255-2 

255-3 

255-5 

255-6 

256-1 

256-2 

256-3 

331 = 

257-2 

257-3 

257-5 

257-6 

258-1 

258-2 

258-3 

258-5 

341 = 

249-2 

249-3 

249-5 

249-6 

250-1 

250-2 

250-3 

250-5 

351= 

243-2 

243-3 

243-5 

243-6 

173-4 

173-5 

173-6 

174-1 

361 = 

174-6 

101-4 

181-5 

181-6 

180-1 

180-2 

180-3 

180-4 

371 = 

179-1 

179-3 

179-4 

179-5 

179-6 

175-1 

175-3 

175-4 

381 = 

177-3 

177-4 

177-5 

177-6 

178-1 

178-2 

178-3 

178-4 

391 = 

176-1 

176-3 

176-4 

176-6 

140-4 

148-5 

148-6 

154 4 


B-4 


- DSP ANALYSIS - 
FREE-FREE MODES 


AUGUST 24, 1903 




ANALYSIS DISPLACEMENT SET 



-1- 

-2- 

-3- 

-4- 

-5- 

-6- 

-7- 

-8- 

1= 

108-1 

108-2 

108-3 

116-1 

116-2 

116-3 

117-1 

117-2 

11 = 

121-3 

132-1 

132-2 

132-3 

123-1 

123-2 

123-3 

106-1 

21 = 

112-1 

112-2 

112-3 

110-1 

110-2 

110-3 

111-1 

111-2 

31 = 

141-1 

141-2 

141-3 

125-1 

125-2 

125-3 

183-2 

185-2 

41 = 

135-3 

143-1 

143-2 

143-3 

139-3 

147-1 

147-2 

147-3 

51 = 

150-3 

269-1 

269-2 

269-3 

272-1 

272-2 

272-3 

255-1 

61 B 

173-2 

173-3 

174-2 

181-1 

181-2 

181-3 

179-2 

175-2 

71 = 

148-1 

140-2 

148-3 

154-1 

154-2 

154-3 

376-1 

376-2 

81 = 

153-2 

153-3 

152-1 

152-2 

152-3 

244-1 

234-1 

228-2 

Sl = 

377-3 

226-2 

156-2 

1-1 

1-2 

1-3 

37-1 

37-2 

101 = 

219-2 

220-1 

220-2 

220-3 

158-2 

218-2 

214-2 

216-2 

111* 

15-2 

15-3 

261-1 

261-2 

261-3 

264-1 

264-2 

264-3 

121 

59-1 

59-2 

59-3 

4-1 

4-2 

4-3 

12-1 

12-2 

131= 

12-5 

12-6 

9-1 

9-2 

9-3 

7-1 

7-2 

7-3 

141 = 

21-3 

28-1 

28-2 

28-3 

97-1 

97-2 

97-3 

73-1 

151= 

199-2 

94-1 

94-2 

94-3 

19-1 

19-2 

19-3 

105-1 

161 = 

40-3 

51-1 

51-2 

51-3 

62-3 

84-1 

84-2 

84-3 

171= 

90-3 

102-1 

102-2 

102-3 

100-1 

100-2 

100-3 

47-1 

181= 

69-1 

69-2 

69-3 

88-1 

88-2 

88-3 

42-3 

44-3 

191 = 

55-2 

55-3 

66-3 







t 


B-5 


- DSP ANALYSIS - 
FREE-i*'REE MODES 


AUGUST 24, 1983 


perm spc Dia;>r , aliment set 


1 = 

-1- 

113-6 

-2- 

1191-5 

-3- 

115-6 

11 = 

189-5 

183-5 

187-5 

21 = 

250-4 

243-4 

174-5 

31 = 

253-4 

244-4 

239-4 

41 = 

226-5 

227-5 

224-5 

51 = 

157-5 

158-5 

162-5 

61 = 

209-5 

210-5 

201-5 

71 = 

50-6 

83-6 

72-6 

81 = 

20-6 

195-5 

197-5 

91 = 

32-6 

93-6 

52-6 

101 = 

77-6 

54-6 



-4- 

-5- 

-6- 

-7- 

193-5 

191-5 

127-6 

133-6 

185-5 

184-5 

144-6 

146-6 

175-5 

176-5 

230-5 

232-5 

165-5 

166-5 

171-5 

169-5 

156-5 

233-4 

245-4 

246-4 

160-5 

218-5 

217-5 

13-6 

242-4 

71-6 

200-5 

204-5 

199-5 

82-6 

237-4 

17-6 

29-6 

B6-6 

74-6 

87-6 

75-6 

92-6 

89-6 

57-6 


1 = 


BDRY SPC DISPLACEMENT SET 
-1- -2- -3- -4- -5- -6- 

255-4 256-4 257-4 234-4 235-4 236-4 


-7- 


► 


« 



l *i 

iS 

3 

$ 

i 


B-G 


-B- 

129-6 

1232-5 

238-4 

167-5 

219-5 

14-6 

208-5 

18-6 

30-6 

80-6 


APPENDIX C 


EXAMPLE OF A STATISTICAL CORRELATION REPORT 


NASTRAN MODAL ANALYSIS - MODAL SURVEY STATISTICAL CORRELATION 
21-N0V-83 

1. INTERACTIVE DIALOG: 

IS A SEPARATE OUTPUT LISTING FILE TO BE PRINTED? (Y OR N) : Y 

ENTER ANALYTICAL LAMA MATRIX FILENAME: LAMA.TBL 
ENTER ANALYTICAL MODE-SHAPE MATRIX FILENAME: PHITE.MTX 
ENTER ANALYTICAL GRID POINT LIST FILENAME! GRDPTF.LIS 

EARNING » UNEQUAL NUMBER OF FREQUENCIES AND MODE SHAPES. 
NUMBER OF FREQUENCIES? 190 

NUMBER OF MODE SHAPES: 48 




ONLY THE 

FIRST 4B 

WILL BE USED. 





ENTER 

NEXT 

EXPERIMENTAL 

MODE- SHAPE 

FILENAME. 

( "NONE" 

IF 

NO 

MORE) 

• 

» 

ENTER 

NEXT 

EXPERIMENTAL 

MODE-SHAPE 

FILENAME. 

( "NONE" 

IF 

NO 

MORE) 

: 

ENTER 

NEXT 

EXPERIMENTAL 

MODE-SHAPE 

FILENAME. 

( "NONE" 

IF 

NO 

MORE) 

• 

• 

ENTER 

NEXT 

EXPERIMENTAL 

MODE-SHAPE 

FILENAME. 

( "NONE" 

IF 

NO 

MORE) 

> 

ENTER 

NEXT 

EXPERIMENTAL 

MODE-SHAPE 

FILENAME. 

( "NONE" 

IF 

NO 

MORE) 

• 

• 

ENTER 

NEXT 

EXPERIMENTAL 

MODE-SHAPE 

FILENAME. 

("NONE" 

IF 

NO 

MORE) 

1 

ENTER 

NEXT 

EXPERIMENTAL 

MODE -SHAPE 

FILENAME. 

( "NONE" 

IF 

NO 

MORE) 

• 

♦ 

ENTER 

NEXT 

EXPERIMENTAL 

MODE-SHAPE 

FILENAME. 

( "NONE" 

IF 

NO 

MORE) 

: 

ENTER 

NEXT 

EXPERIMENTAL 

MODE-SHAPE 

FILENAME. 

("NONE" 

IF 

NO 

MORE) 

s 


IS ANALYTICAL VS. EXPERIMENTAL SYMMETRY TO BE CONSIDERED? <Y OR N) 

PRINT ANALYTICAL MODE-SHAPE VECTORS? (Y OR N) s N 
PRINT EXPERIMENTAL MODE-SHAPE VECTORS? (Y OR N) : N 


RELATIVE DEVIATIONS GREATER THAN A THRESHOLD VALUE WILL BE PRINTED 
THE DEFAULT THRESHOLD IS 0.050 

ENTER DESIRED THRESHOLD. FOR DEFAULT: 0.100 

RELATIVE DEVIATIONS GREATER THAN 10.00 % WILL BE PRINTED. 


C-l 


FMS001.AC2 

FMS002.AC2 

FMS003.AC2 

FMS004.AC2 

FMS005.AC2 

FMS006.AC2 

FMS007.AC2 

FMS008.AC2 

NONE 

N 


ORIGINAL S9f 

OF POOR QUALITY 


NASTRAN MODAL ANALYSIS - MODAL SURVEY STATISTICAL CORRELATION 
21-N0V-83 

2. SUMMARY OF FREQUENCY, MASS, STIFFNESS, DAMPING, AND SYMMETRY : 


ANALYTICAL MODES: 


MODE 

FREQUENCY 

MASS 

STIFFNESS 

SYMMETRY 

1 

4 . 78003426E+01 

1. OOOOOOOOE+OO 

9 . 02031563E+04 

0 

0 

0 

2 

8 . 52370071E+01 

9 . 99999881E-01 

2 , 86824438E+05 

0 

0 

0 

3 

8.92746735E+01 

9 . 99999821E-01 

3 , 14641656E+05 

0 

0 

0 

4 

9 . 368312B4E+01 

9 . 99999940E-01 

3 . 46483469E+05 

0 

0 

0 

5 

1 . 14596474E+02 

1. OOOOOOOOE+OO 

5.18444500E+05 

0 

0 

0 

6 

1 . 29791458E+02 

1.00000000E+00 

6 . 65046375E+05 

0 

0 

0 

7 

1 .42888336E+02 

1 . 00000012E+00 

8.06034125E+05 

0 

0 

0 

8 

1 . 49144562E+02 

1.00000000E+00 

8 , 78161875E+05 

0 

0 

0 

9 

1 . 51106857E+02 

1. OOOOOOOOE+OO 

9, 01422000E+05 

0 

0 

0 

10 

1 . 61486435E+02 

1.00000000E+00 

1.02951300E+06 

0 

0 

0 

11 

1.66986420E+02 

1. OOOOOOOOE+OO 

1 . 10083463E+06 

0 

0 

0 

12 

1.76361176E+02 

1 . 00000012E+00 

1.22790788E+06 

0 

0 

0 

13 

1.86912598E+02 

1. OOOOOOOOE+OO 

1 , 37923075E+06 

0 

0 

0 

14 

2 . 03151978E+02 

9 . 99999940E-01 

1 . 62930288E+06 

0 

0 

0 

15 

2 . 13227921E+02 

1. OOOOOOOOE+OO 

1 , 79493163E+06 

0 

0 

0 

16 

2 . 15393066E+02 

1 . 00000012E+00 

1.83156913E+06 

0 

0 

0 

17 

2 . 20526611E+02 

1 . 00000024E+00 

1 . 91991450E+06 

0 

0 

0 

18 

2 . 25210617E+02 

9 . 99999821E-01 

2 ♦ 0023378QE+06 

0 

0 

0 

19 

2. 28891B00E+02 

1. OOOOOOOOE+OO 

2 . 06833200E+06 

0 

0 

0 

20 

2.37163879E+02 

9 . 99999940E-01 

2 . 22053100E+06 

0 

0 

0 

21 

2 . 41989059E+02 

9 . 99999881E-01 

2 . 31180500E+06 

0 

0 

0 

22 

2 . 44966827E+02 

1. OOOOOOOOE+OO 

2 . 36905050E+06 

0 

0 

0 

23 

2 . 4S436966E+02 

1 , 00000012E+00 

2.39757125E+06 

0 

0 

0 

24 

2 . 59525421E+02 

9.99999881E-01 

2 . 65900725E+06 

0 

0 

0 

25 

2 . 67757904E+02 

1. OOOOOOOOE+OO 

2 . 83037750E+06 

0 

0 

0 

26 

2 . 76532684E+02 

9 . 99999762E-01 

3 . 01892675E+06 

0 

0 

0 

27 

2 . 77095856E+02 

1. OOOOOOOOE+OO 

3 , 03123675E+06 

0 

0 

0 

28 

2.87656219E+02 

1 . 00000012E+00 

3 . 26668575E+06 

0 

0 

0 

29 

2 . 98013702E+02 

1. OOOOOOOOE+OO 

3 , 50616400E+06 

0 

0 

0 

30 

3 . 28465637E+02 

9 . 99999940E-01 

4 . 25931300E+06 

0 

0 

0 

31 

3 . 31687805E+02 

9 . 99999702E-01 

4. 34328750E+06 

0 

0 

0 

32 

3.38717804E+02 

9 . 99999881E-01 

4 . 52934850E+06 

0 

0 

0 

33 

3 . 46994293E+02 

1. OOOOOOOOE+OO 

4 . 75340050E+06 

0 

0 

0 

34 

3 . 55755737E+02 

9 . 99999940E-01 

4 . 99647250E+06 

0 

0 

0 

35 

3 . 83082672E+0'i: 

1 , 00000024E+00 

5.79355150E+06 

0 

0 

0 

36 

3 . 94640564E+02 

1. OOOOOOOOE+OO 

6 . 14841500E+06 

0 

0 

0 

37 

4 . 00549103E+02 

1 , 00000012E+00 

6. 33390250E+06 

0 

0 

0 

38 

4 , O1790157E+O2 

9 . 99999821E-01 

6 . 37346400E+06 

0 

0 

0 

39 

4 . 11575165E+02 

1. OOOOOOOOE+OO 

6 . 68741200E+06 

0 

0 

0 

40 

4 . 23178680E+02 

9 . 99999940E-01 

7 . 06980300E+06 

0 

0 

0 

41 

4 . 62407135E+02 

9 . 99999702E-01 

8 . 44128800E+06 

0 

0 

0 

42 

4 . 66884094E+02 

1. OOOOOOOOE+OO 

8.60553500E+06 

0 

0 

0 

43 

4 . 73970062E+02 

1. OOOOOOOOE+OO 

8 . 86873200E+06 

0 

0 

0 

44 

4 . 76452454E+02 

1 . 00000024E+00 

8 . 96187800E+06 

0 

0 

0 

45 

4. 87654266E+02 

9 . 99999940E-01 

9 , 38823100E+06 

0 

0 

0 

46 

4 . 94913330E+02 

1. OOOOOOOOE+OO 

9 . 66981300E+06 

0 

0 

0 

47 

5 . 04736725E+02 

1. OOOOOOOOE+OO 

1 . 00574890E+07 

0 

0 

0 

48 

5 . 08807892E+02 

1.00000012E+00 

1 . 02203900E+07 

0 

0 

0 
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EXPERIMENTAL MODES i 


MODE FREQUENCY 

DAMPING 

SYMMETRY 

1 

6. 911E+01 

1 • 113E-02 

0 

0 

0 

2 

7 . 811E+01 

4 . 992E-03 

0 

0 

0 

3 

9 . S35E+01 

1 . 076E-02 

0 

0 

0 

4 

I , I30E+02 

1 , 075E-02 

0 

0 

0 

5 

1.645E+02 

4 . 948E-03 

0 

0 

0 

6 

2.258E+02 

3.984E-03 

0 

0 

0 

7 

2 . 280E+02 

2. 466E-03 

0 

0 

0 

8 

3.617E+02 

1.475E-03 

0 

0 

0 
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3. CORRELATION COEFFICIENTS FOR ANALYTICAL VS. EXPERIMENTAL COMPARISONS t 


ANALYTICAL 




EXPERIMENTAL MODES 



MODES 

1 

2 

3 

4 

5 

6 

7 

8 

1 

-0.506 

-0.521 

-0.473 

-0.375 

-0.201 

0.248 

0.327 

0.716 

2 

-0.283 

-0.294 

-0.224 

-0.106 

0.049 

0.324 

0.452 

0.515 

3 

-0.984 

-0.841 

-0.930 

-0.982 

-0.923 

-0.016 

0.016 

0.172 

4 

-0.958 

-0.989 

-0.991 

-0.920 

-0.692 

0.293 

0.338 

0.088 

5 

-0.952 

-0.988 

-0.984 

-0.906 

-0.670 

0.325 

0.372 

0.095 

6 

-0.462 

-0.413 

-0.382 

-0.304 

-0.178 

0.272 

0.404 

0.377 

7 

-0.996 

-0.885 

-0.957 

-0.985 

-0.889 

0.041 

0.078 

0.167 

8 

-0.994 

-0.905 

-0.964 

-0.971 

-0.850 

0.124 

0.166 

0.128 

9 

-0.942 

-0.833 

-0.913 

-0.962 

-0.889 

-0.126 

-0.109 

0.341 

10 

-0.917 

-0.754 

-0.869 

-0.961 

-0.959 

-0.041 

-0.055 

-0.027 

11 

0.971 

0.803 

0.901 

0.966 

0.929 

0.045 

0.011 

-0.197 

12 

0.409 

0.498 

0.523 

0.526 

0.431 

-0.076 

0.007 

0.451 

13 

0.186 

0.408 

0.260 

0.029 

-0.274 

-0.919 

-0.962 

0.091 

14 

-0.035 

-0.121 

-0.149 

-0.181 

-0.176 

0.009 

-0.111 

-0.596 

15 

0.106 

0.186 

0.219 

0.253 

0.238 

0.026 

0.135 

0.575 

16 

0.109 

0.189 

0.222 

0.257 

0.242 

0.063 

0.169 

0.534 

17 

0.722 

0.451 

0.618 

0.803 

0.941 

0.561 

0.545 

-0.173 

18 

0.993 

0.904 

0.965 

0.979 

0.863 

-0.093 

-0.124 

-0.192 

19 

0.961 

0.839 

0.898 

0.910 

0.808 

-0.043 

-0.118 

-0.294 

20 

-0.862 

-0.606 

-0.749 

-0.887 

-0.963 

-0.267 

-0.238 

0.159 

21 

-0.182 

-0.178 

-0.254 

-0.348 

-0.402 

-0.289 

-0.396 

-0.312 

22 

0.078 

0.117 

0.028 

-0.105 

-0.246 

-0.466 

-0.576 

-0.316 

23 

0.227 

0.495 

0.396 

0.231 

-0.038 

-0.715 

-0.647 

0.506 

24 

0.142 

0.198 

0.243 

0.295 

0.297 

0.153 

0.264 

0.364 

25 

-0.135 

-0.167 

-0.224 

-0.296 

-0.324 

-0.224 

-0.334 

-0.337 

26 

-0.248 

-0.389 

-0.275 

-0.088 

0.154 

0.871 

0.912 

-0.037 

27 

0.536 

0.754 

0.668 

0.499 

0.195 

-0.753 

-0.727 

0.341 

28 

-0.109 

-0.137 

-0.197 

-0.272 

-0.310 

-0.215 

-0.327 

-0.400 

29 

-0.212 

-0.490 

-0.343 

-0.109 

0.220 

0.961 

0.958 

-0.326 

30 

-0.214 

-0.471 

-0.324 

-0.092 

0.228 

0.969 

0.979 

-0.257 

31 

-0.230 

-0.516 

-0.359 

-0.111 

0.238 

0.945 

0.964 

-0.255 

32 

0.306 

0.427 

0.431 

0.410 

0.304 

-0.210 

-0.101 

0.397 

33 

0.120 

0.025 

0.008 

-0.029 

-0.043 

-0.108 

-0.306 

-0.262 

34 

0.054 

0.290 

0.206 

0.070 

-0.131 

-0.722 

-0.649 

0.815 

35 

0.248 

0.209 

0.273 

0.356 

0.400 

0.346 

0.431 

-0.336 

36 

-0.426 

-0.420 

-0.426 

-0.416 

-0.333 

-0.152 

-0.149 

0.756 

37 

0.176 

0.409 

0.271 

0.051 

-0.241 

-0,986 

-0.997 

0.359 

38 

-0.916 

-0.958 

-0.939 

-0.841 

-0,592 

0.380 

0.455 

0.070 

39 

0.290 

0.308 

0.339 

0.375 

0.338 

0.473 

0.494 

-0.511 

40 

0.099 

0.029 

0.015 

0.001 

-0.013 

0.270 

0.172 

-0.957 

41 

-0.236 

-0.182 

-0.200 

-0.224 

-0.211 

-0.470 

-0.423 

0.943 

42 

0.133 

0.143 

0.13B 

0.140 

0.096 

0.486 

0.447 

-0.870 

43 

0.153 

0.100 

0.088 

0.075 

0.037 

0.445 

0.318 

-0.950 

44 

0.318 

0.208 

0.265 

0.336 

0.381 

0.371 

0.384 

-0.795 

45 

-0.202 

-0.177 

-0.179 

-0.185 

-0.146 

-0.498 

-0.433 

0.946 

46 

-0.204 

-0.178 

-0.178 

-0.181 

-0.139 

-0.496 

-0.424 

0.957 

47 

0.090 

0.321 

0.193 

0.019 

-0.254 

0.002 

-0.085 

-0.660 

4B 

0.078 

-0.119 

-0.009 

0.128 

0.336 

-0.261 

-0.140 

0.547 
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4. ANALYTICAL MODE SHAPES AND THEIR BEST EXPERIMENTAL MATCHES: 


NASTRAN NASTRAN TEST TEST OF GP'S CORREL MAX REL GRID 

MODE FREQUENCY MODE FREQUENCY COMPARED COEFF DIFFERENC POINT 


1 47.80034 8 361.69601 5 0.716 1.007E+00 234 

RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% i (GRID ID/DEVIATION) 

234-1 /-I . 007E+00 244-1/ 6.232E-01 245-1/ 5.158E-01 255-1/ 8.589E-01 

257-1/ 6 . 611E-01 

SCALED DIFFERENCES ( (XA-XE)/S ) > THRESHOLD i (GRID ID/ DIFFERENCE) 
234-1/-1 . 435E+Q0 244-1/ 8.821E-01 245-1/-7. 546E-01 255-1/ 1.864E+00 

257-1 / -1 . 5B2E+00 


2 85.23701 8 361.69601 5 0.515 1.507E+00 234 

RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% : (GRID ID/DEVIATION) 

234-1 / -1 . 50 7E+00 244-1/-2 . 875E-01 245-1/ 9.712E-01 255-1/ 1.245E+00 

SCALED DIFFERENCES ( <XA-XE)/S ) > THRESHOLD : (GRID ID/ DIFFERENCE) 
234-1/-1. 273E+00 244-1/ 5.324E-01 245-1/-3 . 724E-01 255-1/ 1.536E+00 


j 3 89.27467 1 69.11099 5 -0.984 2.612E-01 257 

} RELATIVE, DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% : (GRID ID/DEVIATION) 

| 244-1/ 1 . 297E-01 255-1 /-2 . 471E-01 257-1 /-2 . 612E-01 

j SCALED DIFFERENCES ( (XA-XE)/S ) > THRESHOLD : (GRID ID/DIFFERENCE) 

! 244-1/ 1 . 253E+02 255-1/ 2.709E+01 257-1/ 3.294E+01 


4 93.68313 3 95.35100 5 -0.991 2.363E-01 234 

RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% : (GRID ID/DEVIATION) 

234-1/ -2. 363E-01 255-1/ 1.033E-01 257-1/ 1.507E-01 


SCALED DIFFERENCES ( (XA-XE)/S ) > THRESHOLD : (GRID ID/DIFFERENCE) 
234-1/ 1 , 430E+01 255-1/ 2.584E+01 257-1/ 2.537E+01 


5 114.59647 2 78.10699 5 -0.988 1.978E-01 234 

RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) ) 10.00% : (GRID ID/DEVIATION) 

234-1 / -1 . 978E-01 244-1/ 1.382E-01 245-1/ 1.853E-01 255-1 /-I . 350E-01 


t ■ 
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SCALED DIFFERENCES ( (XA-XE)/S ) > THRESHOLD i (GRID ID/DIFFERENCE} 
234-1/ 1 . 730E+O1 244-1/ 3.846E+01 245-1/ 3.990E+01 255-1/ 3.214E+C3 


6 129.79146 1 69.11099 5 -0.462 1.871E+00 244 

RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% j (GRID ID/DEVIATION) 

234-1 /-I . 192E+00 244-1/-1. 871E+00 245-1/ 5.432E-01 255-1/-3 . 838E-01 

257-1/ 1 . 195E-01 

SCALED DIFFERENCES ( (XA-XE)/S > > THRESHOLD i (GRID ID/DIFFERENCE) 
234-1 /-3 . 145E+00 244-1/-5. 217E+00 245-1/ 1.041E+01 255-1/ 4.545E-01 

257-1/ 3 . 925E+00 


7 142.8B834 1 69.11099 5 -0.996 1.320E-01 257 

RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% ; (GRID ID/DEVIATION) 

255-1 /-1.081E-01 257-1 /-1.320E-01 

SCALED DIFFERENCES ( (XA-XEJ/S ) > THRESHOLD t (GRID ID/DIFFERENCE) 
255-1/ 6 . 456E+02 257-1/ 7.248E+02 


149.14456 


69.11099 


-0.994 2. 007E-01 


RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% : (GRID ID/DEVIATION) 

234-1/-2. 007E-01 


SCALED DIFFERENCES ( (XA-XEJ/S ) > THRESHOLD : (GRID ID/DIFFERENCE) 
234-1/ 1 . 003E+02 


\ 9 151.10686 4 112.99900 5 -0.962 4.537E-01 255 

i 

l RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% s (GRID ID/DEVIATION) it 

t 234-1/ 1 . 278E-01 245-1/- 1 . 950E-01 255-1/-4 . 537E-01 257-1/ 3.4B1E-01 

£ SCALED DIFFERENCES ( (XA-XE)/S ) > THRESHOLD : (GRID ID/DIFFERENCE) 

234 -1/ -3 . 247E+00 245-1/-4. 100E+00 255-1/-1 . 771E+00 257-1/-1. 181E+00 

■f 

l 10 161.48643 4 112.99900 5 -0.961 3.907E-01 257 

y,, * 

^ RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% i (GRID ID/DEVIATION) 

234-1 / -3 , 256E-01 244-1/ 3.064E-01 255-1 /-1 . 722E-01 257-1/-3 . 907E-01 

SCALED DIFFERENCES ( (XA-XE) /S ) > THRESHOLD : (GRID ID/DIFFERENCE) 

‘ 234-1/-3 . 297E+00 244-1/-4.434E+00 255-1/-1 . 618E+00 257-1/-1 .434E+00 


11 166.98642 1 69.11099 5 0.971 3.647E-01 255 
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RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10,00% i (GRID ID/DEVIATION) 

234-1/ 1 , 162E-01 245-1/-1, 864E-01 255-1/ 3.647E-01 257-1/ 3.218E-01 


SCALED DIFFERENCES ( (XA-XE)/S ) > THRESHOLD s (GRID ID/ DIFFERENCE) 
234-1/-7 . 006E+01 245-1/-1.597E+02 255-1/-1 . 795E+01 257-1/-3. 023E+01 


12 176.36118 4 112.99900 5 0.526 2.016E+00 245 

RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% t (GRID ID/DEVIATION) 

234-1/ 2 . 151E-01 244-1 / -1 . 592E-01 245-1/ 2.016E+00 255-1/-5 . 754E-01 

257-1/ 5 . 240E-01 

SCALED DIFFERENCES ( (XA-XE)/S ) > THRESHOLD : (GRID ID/ DIFFERENCE) 
234-1/ 1 . 046E+00 244-1/ 1.452E+00 245-1/ 1.564E+00 255-1/ 4.019E-01 

257-1/ 4 . 853E-01 


13 186.91260 7 228.02299 5 -0.962 3.708E-01 257 

RELATIVE DEVIATIONS < XA/RMSA-XE/RMSE) > 10.00% j (GRID ID/DEVIATION) 

234-1/-1. 529E-01 244-1/ 1.779E-01 245-1/-2 . 851E-01 255-1/ 3.233E-01 

257-1/-3. 708E-01 

SCALED DIFFERENCES ( (XA-XE)/S ) > THRESHOLD i (GRID ID/DIFFERENCE) 

2 3 4-1/ -4 . 837E+00 244-1/ 6.055E-02 245-1/ 6.729E-01 255-1/ 4.525E+00 

257-1/ 3 . 404E+00 


14 203.15198 8 361.69601 5 -0.596 1.464E+00 234 

RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% t (GRID ID/DEVIATION) 

234-1/ 1 . 464E+00 244-1/ 7.515E-01 245-1/-7 . 9B7E-01 255-1/-5 . 376E-01 

257-1/ 6. 349E-01 

SCALED DIFFERENCES ( <XA-XE)/S > > THRESHOLD : (GRID ID/DIFFERENCE) 
234-1/ 1 . 417E+00 244-1/-2. 671E-01 245-1/ 2.323E-01 255-1/-1 . 368E+00 

257-1/ 1 . 306E+00 


15 213.22792 8 361.69601 5 0.575 1.626E+00 234 

RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% : (GRID ID/DEVIATION) 

234-1/-1. 626E+00 244-1/-7 . 567E-01 245-1/ 6.682E-01 255-1/ 4.587E-01 

257-1 / -6 . 136E-01 


SCALED DIFFERENCES ( <XA-XE)/S ) > THRESHOLD j (GRID ID /DIFFERENCE) 


234-1 / -1 , 616E+00 
257-1 / -1 . 111E+00 

244-1/-1 . 088E-01 

245-1/ 

6 . 512E-02 

255-1/ 1.080E+00 

215.39307 

8 361.69601 

5 

0.534 

1 . 726E+00 234 
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RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% j (GRID ID/DEVIATION) 

234-1/-1 . 7Z6E+00 244-1/-7. 26BE-01 245-1/ 6.917E-01 255-1/ 5.012E-01 

257-1/ -6 . 491E-01 

SCALED DIFFERENCES ( (XA-XE) /3 ) > THRESHOLD i (GRID ID/DIFFERENCE > 
234-1/-1 .431E+00 244-1/ 2.724E-01 245-1/-2.675E-01 255-1/ 1.293E+00 

257-1 / -1 . 250E+00 


17 


220.52661 


164.46001 


0.941 4 , 653E-01 


244 


RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% j (GRID ID/DEVIATION) 

234-1 / -3 . 993E-01 244-1/ 4.653E-01 257-1/ 4.484E-01 


SCALED DIFFERENCES ( (XA-XE) /S ) > THRESHOLD : (GRID ID/DIFFERENCE) 
234-1/-1. 385E+01 244-1/-9. 189E+00 257-1/ 7.887E+00 


18 225.21062 1 69.11099 5 0.993 1.503E-01 244 

RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% : (GRID ID/DEVIATION) 

234-1/ 1 . 279E-01 244-1 / -1 . 503E-01 255-1/ 1.382E-01 


SCALED DIFFERENCES ( (XA-XE) /S ) > THRESHOLD : (GRID ID/DIFFERENCE) 
234-1/ -8 . 470E+01 244-1/-1 . 720E+02 255-1/-5. 054E+01 


19 228.89180 1 69.11099 5 0,961 4.079E-01 245 

RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% j (GRID ID/DEVIATION) 

234-1/ 1.846E-01 244-1/ 3.802E-01 245-1/-4 . 079E-01 255-1/ 2.142E-01 

SCALED DIFFERENCES ( (XA-XE) /S ) > THRESHOLD : (GRID ID/ DIFFERENCE) 
234-1/-3. 699E+01 244-1/-4 . 916E+01 245-1/-1 . 082E+02 255-1/-1 . 950E-* 01 


20 237.16388 5 164.46001 5 -0.963 4.412E-01 245 

RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% : (GRID ID/DEVIATION) 

234-1/-1 . 440E-01 244-1/-3. 758E-01 245-1/ 4.412E-01 


SCALED DIFFERENCES ( (XA-XE) /S ) > THRESHOLD j (GRID ID/DIFFERENCE) 
234-1/ 7 . 447E+00 244-1/ 1.071E+01 245-1/ 1.850E+01 


21 241.98906 5 164.46001 5 -0.402 2.3B3E+00 245 

RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% j (GRID ID/DEVIATION) 

234-1/ 1 , 136E-01 245-1/ -2 . 383E+00 255-1/ 4.441E-01 257-1/-2 . 900E-01 

SCALED DIFFERENCES ( (XA-XE) /S ) > THRESHOLD : (GRID ID /DIFFERENCE) 
234-1/ 1 . 894E+01 245-1/-2 . 310E+01 255-1/ 5.374E+00 257-1/-9 . 868E+00 
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22 244.96603 7 220.02299 5 -0.576 1.290E+00 255 

RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% I (GRID ID/DEVIATION) 

234-1/-4. 072E-01 244-1/ 1.069E+00 245-1/-1 . 123E+00 255-1/ 1.290E+00 

SCALED DIFFERENCES ( <XA-XE)/S ) > THRESHOLD i (GRID ID/DIFFERENCE) 
234-1/-1 .553E+00 244-1/ 2.123E-01 245-1/ 1.399E-02 255-1/ 1.615E+00 


23 246.43697 6 225.84001 5 -0.715 1.142E+00 234 

RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% i (GRID ID/DEVIATION) 

234 -1/ -1 . 142E+00 244-1/-8.476E-01 245-1/ 8.765E-01 255-1/-2. 123E-01 

257-1/ 1.105E-01 

SCALED DIFFERENCES ( (XA-XE)/S ) > THRESHOLD i (GRID ID/DIFFERENCE) 
234-1/ -2 . 057E+00 244-1/-1. 216E-01 245-1/ 3.905E-01 255-1/ 1.124E+00 

257-1/ 1 . 002E+00 


24 259.52542 0 361.69601 5 0.364 1.095E+OO 234 

RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% i (GRID ID/DEVIATION) 

234-1/-1 . 895E+00 244-1/-8.343E-01 245-1/ 7.344E-01 255-1/ 0.61OE-O1 

257-1/ -8 . 807E-O1 

SCALED DIFFERENCES ( <XA~XE)/S ) > THRESHOLD j (GRID ID/DIFFERENCE) 
234-1/-1 . 635E+00 244-1/-2 .462E-01 245-1/ 1.926E-01 255-1/ 1.135E+00 

257-1/-1 . 103E+00 


25 267.75790 8 361.69601 5 -0.337 1.970E+00 234 

RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% i (GRID ID/DEVIATION) 

234-1/ 1 . 970E+00 244-1/ 8.220E-01 245-1/-7. 023E-01 255-1/-9 . 740E-01 

257-1/ 7 . 927E-01 

SCALED DIFFERENCES ( (XA-XE)/S ) > THRESHOLD t (GRID ID /DIFFERENCE) 
234-1/ 2 . 7B1E+00 244-1/ 1.894E+00 245-1/-1 . 668E+00 255-1/-7 . B08E-01 

257-1/ 5 . 379E-01 


26 276.53268 7 220.02299 5 0.912 6.383E-01 245 

RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% s (GRID ID/DEVIATION) 

245-1/ 6 . 303E-O1 255-1/-6. 327E-01 257-1/ 2.610E-01 

SCALED DIFFERENCES ( (XA-XE) /S ) > THRESHOLD i (GRID ID/ DIFFERENCE) 
245-1/-4. 325E-01 255-1/-3 . 131E+00 257-1/-2 . 387E+00 


27 277.09586 2 78.10699 5 0.754 1.107E+00 234 
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RELATIVE DEVIATIONS ( XA/RMSA-XE/RMSE) > 10.00% t (GRID ID/DEVIATION) 

234-1/ 1. 107E+00 244-1/-1 . 991E-Q1 245-1/ 1.035E+00 255-1/-3.451E-01 

SCALED DIFFERENCES ( (XA-XE)/S ) > THRESHOLD t (GRID ID/DIFFERENCE) 
234-1/ 4 . 357E+00 244-1/-9. 549E+00 245-1/ 1.9B1E+00 255-1/-1 . 196E+01 


28 287.65622 8 361.69601 5 


-0,400 1, 903E+00 234 


RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% r (GRID ID/DEVIATION) 

234-1/ 1 , 903E+00 244-1/ 8.008E-01 245-1/-6. 998E-01 255-1/-9 . 072E-01 

257-1/ 6. 519E-01 


SCALED DIFFERENCES ( (XA-XE)/S ) > THRESHOLD » (GRID ID/DIFFERENCE) 
234-1/ 4 . 232E+00 244-1/ 4.260E+00 245-1/-3 . 861E+00 255-1/ 2.393E-02 

257-1/-6.488E-01 


29 298.01370 6 225.84801 


0.961 4 . 744E-01 234 


RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% : (GRID ID/DEVIATION) 

234-1/ 4 , 744E-01 244-1/ 1.119E-01 255-1/ 1.538E-01 257-1/ 3.453E-01 


SCALED DT'FFERENCES ( (XA-XE)/S ) > THRESHOLD i (GRID ID/DIFFERENCE) 
234-1/ 1 . 703E+00 244-1/ 4.127E-01 255-1/ 5.829E-01 257-1/ 1.277E+00 


30 328.46564 7 ?.?. 8. 02299 5 0.979 3.031E-01 257 

RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% : (GRID ID/DEVIATION) 

234-1/ 2 . 192E-01 244-1/ 2.321E-01 257-1/ 3.031E-01 


SCALED DIFFERENCES { (XA-XE)/S ) > THRESHOLD t (GRID ID/DIFFERENCE) 


234-1/ 2 . 075E+00 

244-1/ 9.637E-01 

257-1/ 3 , 996E-01 


331.68781 

7 228.02299 

5 0.964 4 . 343E-01 

234 


RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% : (GRID ID/DEVIATION) 

234-1/ 4.343E-01 244-1/ 1.634E-01 245-1/-1. 717E-01 255-1/ 1.0UE-01 

257-1/ 3 . 285E-01 


SCALED DIFFERENCES ( (XA-XEJ/S ) > THRESHOLD t (GRID ID/DIFFERENCE) 
234-1/ 4.766E+00 244-1/ 1.005E-01 245-1/-8 , 296E-01 255-1/-4 . 088E+00 

257 -1/ -3 . 117E+00 


32 338.71780 3 95.35100 5 0.431 2.228E+00 245 

RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% j (GRID ID/DEVIATION) 

234-1/ 4 . 634E-01 244-1/-4 . 797E-01 245-1/ 2.228E+00 255-1/ 1.077E-01 

257-1/ 5. 222E-01 
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SCALED DIFFERENCES ( (XA-XE) /S ) > THRESHOLD i {GRID ID/DIFFERENCE) 
234-1/ -5. 728E-02 244-1/-2 . 912E+00 245-1/ 4.001E+00 255-1 /-9 . 730E-01 

257-1/ 1.745E-01 


33 346.9942 9 7 228.02299 5 -0.306 1.561E+00 245 

i * RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% « (GRID ID/DEVIATION) 

| 234-1/-6. 545E-01 244-1/ 1.064E+00 245-1/-1 . 561E+00 255-1/ 1.459E+00 

J 257-1/ 9 , 013E-01 

' SCALED DIFFERENCES ( (XA-XE) /S ) > THRESHOLD « (GRID ID/ DIFFERENCE) 

234-1 / -1 . 449E+00 244-1/ 8.250E-O2 245-1/ 1.038E-01 255-1/ 1.403E+00 

257-1/ 1 . 101E+00 


34 355.75574 8 361.69601 5 0.815 1.242E+00 257 

RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% : (GRID ID/DEVIATION) 

234-1/ -2 . 168E-01 244-1/ -2 . 589E-01 245-1/ 2.124E-01 255-1/-3.909E-01 

257-1 / -1 , 242E+00 

SCALED DIFFERENCES ( (XA-XE) /S ) > THRESHOLD : (GR.D ID/ DIFFERENCE) 
234-1 /-I . 312E+0Q 244-1/ 5.981E-01 245-1/-5 . 989E-01 255-1/ 1.383E+00 

257-1/ -2 . 116E+00 


35 383.08267 7 228.02299 5 0.431 1.747E+00 257 

RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% i (GRID ID/DEVIATION) 

234-1/ 3 . 683E-01 244-1/-1 . 352E+00 245-1/ 6.752E-01 255-1/-4. 676E-01 

257-1 / -1 , 747E+00 

SCALED DIFFERENCES ( (XA-XE' '*! ) > THRESHOLD i (GRID ID/DIFFERENCE) 
234-1/ 1 . 523E+00 244-1/ -8 . 820E-02 245-1/-1 . 971E-01 255-1/ -1 . 413E+00 

257-1 /-I . 220E+00 


36 394 .6-t056 8 361.69601 5 0.756 1.308E+00 244 

RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% : (GRID ID/DEVIATION) 

234-1/-1. 235E-01 244-1/ 1.308E+00 245-1/-5. 998E-01 255-1/ 3.757E-01 

257-1/ 4. 632E-01 

SCALED DIFFERENCES < (XA-XE) /S ) > THRESHOLD s (GRID ID /DIFFERENCE) 
234-1 /-l . 308E+00 244-1/ 1.115E+00 245-1/-9 . 107E-01 255-1/ 1.792E+00 

257-1 1„4' , 9E+00 


37 54910 7 228.02299 5 -0.997 1.568E-01 244 

RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% t (GRID ID/DEVIATION) 
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244-1/-1. 5681-01 

SCALED DIFF ERENCES ( (XA-XE) /S } > THRESHOLD t (GRID ID/DIFFERENCE) 
244-1/ -1.598E -01 


38 401.79816 2 78.10&W r j -0.958 4.568E-Q1 245 

RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% s (GRID ID/ DEVIATION) 

2 34-1/ -3 . 745E-01 244-1/-2 . 143E-01 545-1/ 4.568E-01 257-1/-1 . 688E-01 

SCALED DIFFERENCES ( (XA-XE) /S ) > THRESHOLD j (Gf.ID ID/DIFFERENCE) 
234-1/-2.346E-01 244-1/ 1.404E+00 245-1/ 5.543E+00 257-3/ 1. 934,3+1/0 


39 411.57516 8 361.69601 % -0.511 1.362E+00 245 

RELATIVE DEVIATIONS < XA/RMSA-XE/RMSE) > 10.00% : (GRID ID/DEVIATION' 

234-1/ -6 , 936E-01 244-1/-1.345E+00 245-1/ 1.362E+00 255-1/-7. 729E-01 

257-1/ 3 . 877E-01 

SCALED DIFFERENCES ( (XA-XE) /S ) > THRESHOLD : (GRID ID/ DIFFERENCE) 
234-1/ 4 . 946E-01 244-1 / -1 . 005E+00 245-1/ 9.959E-01 255-1/ -1 . 32BE+00 

257-1/ 1 , 090E+00 


40 423.17868 8 361.69601 5 -0.957 5.475E-01 234 

RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% : (GRID ID/DEVIATION) 

234-1/ 5.475E-01 244-1/ 2.781E-01 245-1/-2 . 173E-01 

SCALED DIFFERENCES ( (XA-XE) /S ) > THRESHOLD : (GRID ID /DIFFERENCE) 
234-1/ 2.881E+00 244-1/-9 . 786E-01 245-1/ 1.014E+00 


43. 462.40714 8 361.69601 5 0.943 5.118K~01 244 

RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% s (GR.D ID/DEVIATION) 

234-1/ 2 . 562E-01 244-1/ 5.118E-01 245-1/-4 . 861E-01 

SCALED DIFFERENCES ( (XA-XE) /S ) > THRESHOLD : (GRID ID/DIFFERENCE) 
234-1/ 2 . 506E+00 244-1/ 1.391E+00 245-1/-1 .306E+00 


42 466.88409 8 361.69601 5 -0.870 9.348E-01 245 

RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% : (GRID ID/DEVIATION) 

234-1/-4. 404E-01 244-1/-4. 085E-01 245-1/ 5.348E-01 255-1/-2. 553E-01 

SCALED DIFFERENCES ( (XA-XE) /S ) > THRESHOLD i (GRID ID/DIFFERENCE) 
234-1/-9.477E-01 244-1/-8. 201E-01 245-1/ 1.908E+00 255-1/-4 . 763E-01 
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43 473.97006 8 361.69601 5 -0.950 5.011E-01 244 


RELATIVE DEVIATIONS ( J^A / RMP-A-XE / RMSE ) > 10.00% i (GRID ID/DEVIATION) 

244-1/ 5 . 011E-01 245-1/-2 .'25QE-01 255-1/-3 . 669B-01 257-1/ 2.301E-01 

SCALED DIFFERENCES ( (XA-XE)/S ) > THRESHOLD : (GRID ID/DIFFERENCE) 
244-1/ -1 . 464E-02 245-1/ 4.715E-01 255-1/-2 . 585E+00 257-1/ 2.183E+00 


44 476.45245 8 361.&9601 5 -0.795 1.302E+00 244 

RELATIVE DEVIATIONS {XA/RMSA-XE/RMSE) > 10.00% : (GRID ID/DEVIATION) 

234-1/ 1 . 234E-01 244-1/-1.302E+00 245-1/ 5.608E-01 255-1/ 1.600E-01 

SCALED DIFFERENCES ( (XA-XE)/S ) > THRESHOLD : (GRID ID/ DIFFERENCE) 
234-1/ 1.174E+00 244-1/-1 . 387E+00 245-1/ 9.626E-01 255-1/-1 . 384E+00 


45 487.65427 0 361.69601 5 0.946 5.167E-01 245 

RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% s (GRID ID/DEVIATION) 

234-1/ 4.123E-01 244-1/ 2.499E-01 245 ,/~5. 167E-01 255-1/ 2.058E-01 

SCALED DIFFERENCES ( (XA-XE)/S ) > THRESHOLD i (GRID ID/DIFFERENCE) 
234-1/ 2 . 616E-f-00 244-1/ 4.199E-01 245-1/-1 . 514E+00 255-1/-4.404E-01 

V 


46 494.91333 8 361.69601 5 0.957 4.386E-01 245 

RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% i (GRID ID/DEVIATION) 

234-1/ 4.023E-01 244-1/ 1.726E-01 245-1/-4 . 3S6E-01 255-1/ 2.224E-01 

SCALED DIFFERENCES ( (XA-XE)/S ) > THRESHOLD : (GRID ID/DIFFERENCE) 
234-1/ 1 . 448E+01 244-1/ -3 . 542E+00 245-1 /-3 . 938E-02 255-1/-9.586E+00 


47 504.73672 8 ,361.69601 5 -0.660 1.438E+00 255 

RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% : (GRID ID/DEVIATION) 

244-1/ 1 . 999E-01 245-1/ 1.002E+00 255-1/-1 . 43BE+00 257-1/-5. 369E-01 

SCALED DIFFERENCES ( (XA-XE)ZS > > THRESHOLD t (GRID ID/DIFFERENCE) 
244-1/-3 . 375E-01 245-1/ 1.033E+00 255-1/-1 . 835E+00 257-1/ 5.864E-01 


48 508.80789 8 361.69601 5 0.547 1.649E+00 255 

RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% : (GRID ID/DEVIATION) 

234-1/ 3 . 762E-01 244-1/-8 . 422E-01 245-1/-9 . 169E-01 255-1/ 1.649E+00 

257-1 / -3 . 427E-01 
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SCALED DIFFERENCES ( (XA-XE)/S > > THRESHOLD : (GRID ID/DIFFERENCE) 
234-1/ -8 , 649E-0X 244-1/ 3.842E-01 245-1/-7. 889E-01 255-1/ 1.667E+00 

257-1/-1 . 247E+00 


NASTRAN MODAL ANALYSIS - MODAL SURVEY STATISTICAL CORRELATION 
21-N0V-83 

5. EXPERIMENTAL MODE SHAPES AND THEIR BEST ANALYTICAL MATCHES: 

TEST TEST NASTRAN NASTRAN OF GP'S CORREL MAX REL GRID 

MODE FREQUENCY MODE FREQUENCY COMPARED COEFF DIFFERENC POINT 


1 69.11099 7 142.88834 5 -0.996 1.320E-01 257 

RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% : (GRID ID/DEVIATION) 

255-1/-1 . 081E-01 257-1 /-I . 320E-01 

SCALED DIFFERENCES ( (XA-XE)/S ) > THRESHOLD : (GRID ID/DIFFERENCE) 
255-1/ 6 . 456E+02 257-1/ 7.248E+02 


2 78.10699 4 93.68313 


-0.989 1 . 799E-01 244 


RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% : (GRID ID/DEVIATION) 

2 34-1 /-I . 252E-01 244-1/ 1.799E-01 245-1/ 1.545E-01 255-1 /-I . 372E-01 

257-1 / -1 . 226E-01 


3CALED DIFFERENCES < <XA-XE)/S ) > THRESHOLD : (GRID ID/DIFFERENCE) 
234-1/ 3.614E+01 244-1/ 6.954E+01 245-1/ 6.712E+01 255-1/ 5.687E+01 

257-1/ 5 . 509E+01 


3 95.35100 4 93.68313 5 -0.991 2.363E-01 234 

RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% : (GRID ID/DEVIATION) 

234-1/-2. 363E-01 255-1/ 1.033E-01 257-1/ 1.507E-01 


SCALED DIFFERENCES ( (XA-XE)/S ) > THRESHOLD r (GRID ID/DIFFERENCE) 
234-1/ 1 . 430E+01 255-1/ 2.584E+01 257-1/ 2.537E+01 


4 112.99900 7 142.B8B34 5 -0.985 2.519E-01 245 

RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% : (GRID ID/DEVIATION) 

234-1/-2. 369E-01 244-1/-1 . 197E-01 245-1/ 2.519E-01 257-1/ 1.145E-01 

SCALED DIFFERENCES ( (XA-XE)/S ) > THRESHOLD : (GRID ID/DIFFERENCE) 
234-1/ 2 . 785E+00 244-1/ 6.626E+00 245-1/ 9.541E+00 257-1/ 3.-I60E+00 


5 164.46001 20 237.16388 5 -0.963 4.412E-01 245 

RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% : (GRID ID/DEVIV’ r! ‘GN> 

234-1/-1. 440E-01 244-1/-3 . 758E-01 245-1/ 4.412E-01 
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SCALED DIFFERENCES < (XA-XE)/S > > THRESHOLD i (GRID ID/DIFFERENCE) 
234-1/ 7 . 447E+00 244-1/ 1.071S+01 245-1/ 1.850E+01 


6 225.84001 37 400.54910 5 -0.986 2.857E-01 245 

RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% s (GRID ID/DEVIATION) 

234-1/-1 . 719E-01 244-1/ 1.032E-01 245-1/-2 . 857E-01 257-1/-1 . 082E-01 

SCALED DIFFERENCES ( (XA-XE)/S ) > THRESHOLD i (GRID ID /DIFFERENCE) 
234-1 / -8 . 159E+00 244-1/ 1.411E+00 245-1/-2 . 381E-01 257-1/ 4.829E+00 


7 228.02299 37 400.54910 5 -0.997 1.568E-01 244 

RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% i (GRID ID/DEVIATION) 

244-1/-1. 568E-01 

SCALED DIFFERENCES ( (XA-XEJ/S > > THRESHOLD j (GRID ID/DIFFERENCE) 
244-1/-1. 598E-01 


8 361.69601 40 423.17868 5 -0.957 5.475E-01 234 

RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > 10.00% : (GRID ID/DEVIATION) 

234-1/ 5. 475E-01 244-1/ 2.781E-01 245-1/-2. 173E-01 

SCALED DIFFERENCES ( (XA-XEJ/S ) > THRESHOLD s (GRID ID /DIFFERENCE) 
234-1/ 2 . 881E+00 244-1/-9. 786E-01 245-1/ 1.014E+00 


APPENDIX D 


FORTRAN SOURCE CODE LISTING (TESETDMI .FOR) 


PROGRAM TESETDMI 
C 

C THIS PROGRAM WILL GENERATE 'THE TESET MATRIX DMI INPUT TO INSERT INTO 
C THE NASTRAN BULK DATA, FOR USE IN PARTITIONING THE PHIG MATRIX TO 
C OBTAIN THE PHITE MATRIX. 

C INPUT FILES : 

C EXTERNAL GRIDPOINT AND DOF FILE - CONTAINS THE GRIDPOINT ID's 

C AND THEIR ASSOCIATED DOF COMPONENTS. CREATED BY THE USER WITH 

C ALL ENTRIES IN LIST-DIRECTED (FREE-FIELD) FORMAT. 

C ENTRIES ARE INSERTED IN GPID,DOF PAIRS, ONE PAIR FOR EACH 

C GRIDPOINT AND DOF-COMPONENT WHICH IS INSTRUMENTED IN THE 

C MODAL SURVEY. 

C 

C NASTRAN PRINT FILE WITH DTAG 21 OUTPUT. 

C 

PARAMETER < MAXMSDOF* 10 00 > I MAXIMUM NO. OF MODAL SURVEY GPID-DOF's 
C 

INTEGER GPLIST( MAXMSDOF), GPDOF ( MAXMSDOF > , INTDOF ( MAXMSDOF ) 

INTEGER IDOF, EGP , EDOF, G 
CHARACTER REC*120, MSG211BA52, ENDOF21*37 
CHARACTER' 1 ' 40 D21FILE, MSGPFILE, DMIFILE 
C 

DATA M5G2118 /'*** USER INFORMATION MESSAGE 2118, SUBROUTINE GP4PRT' / 

DATA ENDOF21 /' COLUMN TOTALS '/ 

C 

C GET NASTRAN DIAG 21 PRINT FILE 
TYPE 1000, 'FILENAME OF NASTRAN PRINT FILE WITH DIAG 21 i' 

ACCEPT 2000, D21FILE 

OPEN { UNIT- 2 , NAME=D2 1FILE ,TYPE= ' OLD ' , READONLY ) 

C 

C GET MODAL SURVEY GRID-POINT LIST FILE 
TYPE 1000, 'FILENAME OF MODAL SURVEY GRIDPOINT LIST:' 

ACCEPT 2000, MSGPFILE 

OPEN ( UNIT* 1 , NAME=MSGPFILE , TYPE" ' OLD ' , READONLY ) 

C 

C GET TESET DMI FILENAME 

TYPE 1000, 'FILENAME FOR TESET DMI BULK DATA RECORDS:' 

ACCEPT 2000, DMIFILE 
C 

C READ MODAL SURVEY GPID-DOF LIST 
READ(1,*,END=10) ( GPL 1ST ( I > ,GPDOF( I ), 1=1 , MAXMSDOF) 

10 CLOSE(UNIT=l ) 

C 

C COUNT AND PACK MODAL SURVEY GPID-DOF'S 
1 = 1 

DO WHILE (GPLIST(I) .NE. 0) 

GPLIST(I) = lOAGPLIST(I) + GPDOF(I) 

I = 1+1 

IF (I .GT. MAXMSDOF) GOTO 20 


D-l 


NDOFS « 1-1 


END DO 
20 
C 

C SORT GPID-DOF LIST 
DO 100 I°1 ,NDOFS-l 
DO 100 J«I+1, NDOFS 

IF (GPLIST(J) .LT. GPLIST(D) 

1 CALL SWAP(GPLIST{ I) , GPLIST(J)) 

100 END DO 

C 

C FIND DIAG 21 IN PRINT FILE 
DO WHILE (REC( 1 3 52 ) .NE. MSG2118) 

READ( 2, 3000, END* 910) REC 
END DO 
C 

C MATCH INTERNAL DOF'S (FROM DIAG 21) WITH MODAL SURVEY GPID-DOF'S 
DO 300 WHILE <REC(1j37) .NE. ENDOF21 ) 

READ(2,3000,END=920) REC 
IF (REC( 35 i 35 ) .NE. '-') GOTO 300 
READ(REC,4000 ,ERR a 300 ) IDOF, EGP, EDOF 
G *« IDOF 

EGP = 10*EGP + EDOF IPACK G.P. AND DOF 
C CHECK FOR EGP IN MODAL SURVEY LIST 

CALL BINSRCH( EGP , GPLIST , NDOFS , LOC ) 

IF (LOC .NE. 0) INTDOF(LOC) a IDOF 
300 END DO 

CL0SE(UNIT=2) 

C 

C SORT LISTS ON INTERNAL DOF VALUE 
DO 400 I -1, NDOFS -l 
DO 400 Jel+i, NDOFS 

IF (INTDOF(J) .LT. INTDOF(I)) THEN 
CALL SWAP( INTDOF(I), INTDOF(J) ) 

CALL SWAP( GPLIST ( I) , GPLIST(J) ) 

END IF 

400 END DO 

C 

C UNPACK GPID-DOF'S 
DO N=l, NDOFS 

GPDOF(N) a MOD( GPLIST(N), 10 > 

GPLIST(N) = GPLIST(N) / 10 
END DO 
C 

C CHECK FOR G.P.'S WITHOUT INTERNAL DOF MATCH 
N = 1 

DO WHILE ( ( INTDOF(N) .EQ. 0) .AND. (N .LE. NDOFS)) 

TYPE 5000, GPLIST! N) , GPDOF(N) 

N = N+l 

IF (N .GT. MAXMSDOF) GOTO 500 
END DO 

500 IF (N .GT. NDOFS) THEN 

TYPE A, 'ERROR - NO INTERNAL DOF FOR ANY MODAL SURVEY POINT' 


STOP 'FORTRAN STOP — EXECUTION ABORTED' 

END IF 
N START = N 
C 

C WRITE DMI BULK DATA RECORDS 

OPEN(UNIT=3,NAME=DMIFILE,TYPE='NEW' , CARRl'AGECONTROL= ' LIST ' ) 
WRITE(3,6000) 0,2,1,1,G,1 
ISEQ - 1 

IF ( NDOFS .GT. NSTART+2 ) THEN 
WRITE( 3 ,7000 ) 1, 

1 (INTDOF(N) ,GPLIST(N) ,GPDOF(N) ,N=NSTART,NSTART+2 ) , ISEQ 

ELSE 

NT = NDOFS - NSTART + 1 

WRITE( 3 , 8000 ) 1,(INTD0F(N) ,GPLIST(N) ,GPDOF(N) ,N=NSTART,NDOF) 

END IF 

I = NSTART + 3 
DO WHILE (I ,LE. NDOFS) 

IF (NDOFS .GT. 1+3) THEN 
WRITE( 3 ,9000 ) ISEQ, 

1 (INTDOF(N) ,GPLIST(N) ,GPDOF(N) ,N=I,I+3) , ISEQ+1 

ISEQ = ISEQ + 1 

ELSE 

NT = NDOFS - I + I 
WRITE( 3 , 10000 ) ISEQ, 

I (INTDOF(N) ,GPLIST(N) ,GPDOF(N> ,N=I, NDOFS) 

END IF 
1 = 1 + 4 
END DO 

CL0SE(UNIT=3 ) 

C 

STOP 'FORTRAN STOP — PROCESSING COMPLETED' 

C 

910 TYPE *, 'NO DIAG 21 OUTPUT FOUND IN NASTRAN PRINT FILE.' 

STOP 'FORTRAN STOP — EXECUTION ABORTED' 

C 

920 TYPE *, 'END-OF-FILE WHILE READING NASTRAN DIAG 21 OUTPUT' 

STOP 'FORTRAN STOP — EXECUTION ABORTED' 

C 

1000 FORMAT ( IX, A, ' ',$) 

2000 FORMAT ( A ) 

3000 FORMAT ( IX ,A120) 

4000 FORMAT ( 12 2, 111, 2X, 12) 

5000 FORMAT ( ' WARNING - NO INTERNAL DOF FOR GRID POINT' ,19, ,11) 

6000 FORMAT ( 'DMI ' ,5X, 'TESET' ,3X,4I8,8X,2I8) 

7000 FORMAT ( 'DMI ' ,5X, 'TESET' ,3X,I8,3 ( 18 ,16, ' . ' ,11) , '+TE' ,15.5) 

8000 FORMAT ( ' DMI ' ,5X, 'TESET' ,3X,IB, <NT>( 18,16, ' . ' ,11) ) 

9000 FORMAT ( ' +TE' ,15.5,4(18,16, ' . ' ,11) ,'+TE' ,15.5) 

10000 FORMAT ( '+TE' , 15 . 5 , <NT> ( 18,16, ' . ' ,11)) 

END 
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SUBROUTINE BINSRCH ( IV, LIST ,L,LOC) 

C 

C PERFORM BINARY SEARCH FOR VALUE IV IN LIST. 
C LIST MUST BE IN ASCENDING SORT. 

C LOCATION OF IV IS RETURNED IN LOC. 

C IF IV IS NOT FOUND, L0C=0 IS RETURNED. 

C 

DIMENSION LIST(*) 

C 

LOC = U 

IF <L .LT. 1) RETURN 
LBTM - 1 
LTOP = L 
C 

DO WHILE (LTOP-LBTM .GT. 1) 

LMID = (LBTM+LTOP) / 2 
IF (IV .LT. LIST(LMID) ) THEN 
LTOP = LMID 

ELSE IF (IV .GT. LIST(LMID)) THEN 
LBTM = LMID 

ELSE 

LOC = LMID 
RETURN 
END IF 
END DO 
C 

IF (IV .EQ. LIST (LBTM) ) THEN 
LOC = LBTM 

ELSE IF (IV .EQ. LIST(LTOP)) THEN 
LOC = LTOP 
END IF 
C 

RETURN 

END 


SUBROUTINE SWAP(I.J) 
C 

K=I 

I=J 


J=K 


RETURN 

END 


APPENDIX E. 


FORTRAN SOURCE CODE LISTING (LAMA. FOR) 


PROGRAM LAMATABLE 
C 

PARAMETER (MAXEIG = 300) 

C 

CHARACTER LAMA^B 

REAL FREQ ( MAXEIG ) , GMASS ( MAXEIG ) , STIFF ( MAXEIG ) , S YMM ( 3 ) 

REAL REC2(7*MAXEIG) 

C 

DATA LUNP, LUNM fl , 2/ 11 PUNCH AND MATRIX FILE LOG UNIT 'S 

C 

C POSITION PUNCH FILE TO 1ST CARD OF RECORD 2 OF LAMA TABLE 
CALL FINDREC2( LUNP, LAMA) 

C READ RECORD 2 OF LAMA TABLE INTO ARRAY REC2 
CALL READREC2 ( REC2 ,NVALUES , LUNP , MAXEIG ) 

C EXTRACT FREQUENCY, MASS, AND STIFFNESS VECTORS FROM REC2 
NEIGENS = NVALUES / 7 
DO 100 1=1, NEIGENS 

FREQ(I) = REC2(7*I~2) 

GMASS ( I ) = REC2 ( 7*1-1 ) 

STIFF(I) = REC2(7*I) 

100 CONTINUE 

C DO ASCENDING SORT OF FREQ VECTOR? SLAVE SORT MASS AND STIFF VECTORS 
CALL SORT ( NEIGENS , FREQ , GMAS S , STIFF ) 

C CHECK PUNCH FILE FOR SYMMETRY PARAMETERS 

CALL READS YM(S YMM, LUNP) 

C CLOSE PUNCH FILE 

CLOSE ( UNIT=LUNP ) 

C WRITE MATRIX FILE 

CALL WRITEMTX ( NEIGENS , FREQ , GMASS , STIFF , S YMM , LUNM , LAMA ) 

STOP 

END 
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SUBROUTINE FINDREC2 ( LUNP , LAMA ) 

C 

CHARACTER PFILEA80 ,LAMA*B,LINE*80 
C 

C OPEN PUNCH FILE 

10 TYPE 1000, 'ENTER INPUT PUNCH FILE NAME;' 

ACCEPT 2000, PFILE 

OPEN ( UNIT=LUNP , NAME=PFILE , TYPE* ' OLD ' , READONLY , ERR* 1 0 ) 

C POSITION FILE TO LAMA TABLE 

TYPE 1000, 'ENTER NAME OF LAMA TABLE j ' 

ACCEPT 3000, LAMA 
DO 100 1=1,1000000 

READ (LUNP, 200 0,END=1 25) LINE 

IF ( (LINE(ls4> .EQ. 'DTI*' ) .AND. (LINE( 9 j16) .EQ.LAMA) ) 

A GOTO 150 II LAMA TABLE FOUND 

100 CONTINUE 

C TABLE NOT FOUND; START OVER 

125 TYPE A, 'A AAA LAMA TABLE NOT FOUND' 

CLOSE (UNIT=LUNP) 

GOTO 10 

C POSITION PUNCH FILE TO RECORD 2 OF LAMA TABLE 
150 DO 200 1=1,1000000 

READ( LUNP, 4000 ,END=225) LINE<1;24) ,IREC 

IF ( (LINE(1:4) .EQ. 'DTIA' ) .AND. (LINE! 9 i 16) .EQ.LAMA) 

A .AND. (IREC.EQ.2)) GOTO 250 11 RECORD 2 FOUND 

200 CONTINUE 

C RECORD 2 NOT FOUND; START OVER 

225 TYPE A, ' AAA* RECORD 2 OF LAMA TABLE NOT FOUND' 

CLOSE! UN IT=LUNP) 

GOTO 10 

C POSITION FILE TO REREAD 1ST CARD OF RECORD 2 

250 BACKSPACE LUNP 

RETURN 
C 

1000 FORMAT! IX, A30, IX, $> 

2000 FORMAT (A80) 

3000 FORMAT! A8) 

4000 FORMAT! A24, 116) 

END 
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SUBROUTINE READREC2 ( I REC2 , NVALUES , LUNP , MAXEIG ) 

C 

INTEGER IREC2(7*MAXEIG),ITEMP(4) 

CHARACTERS CTEMP(4) 

C 

C READ 1ST CARD OF LAMA TABLE RECORD 2 
READ (LUNP ,1000) IREC2( 1 ) , IREC2 ( 2) 

NVALUES = 2 

C READ REMAINING CARDS UNTIL 'ENDREC' FOUND 
DO 200 1=1,1000000 

READ( LUNP, 2000 ,END= 225 ) (CTEMP< J) ,ITEMP( J) ,J=1,4) 

DO 100 J=1 ,4 

IF (CTEMP(J) .EQ. 'ENDREC') GOTO 250 
NVALUES = NVALUES + 1 
IREC2< NVALUES) = ITEMP(J) 

IF (CTEMP(J) (6:6) .EQ. 

IREC2( NVALUES) = -IREC2 (NVALUES ) 

IF (NVALUES .EQ. 7 >*< MAXEIG) GOTO 225 
CONTINUE 
CONTINUE 

END OF RECORD 2 OF LAMA TABLE NOT FOUND 
TYPE 3000, NVALUES 
PAUSE 'TYPE "CONTINUE" OR "STOP"' 

0 CHECK TO SEE IF OF RECORD 2 ENTRIES IS DIVISIBLE BY 7 
250 IF (MOD (NVALUES, 7) .NE. 0) TYPE 4000, NVALUES 

RETURN 
C 

1000 FORMAT(40X,2I16) 

2000 FORMAT ( 8X , 4 ( A6 ,110 ) ) 

3000 FORMAT ( IX, ' A** WARNING j END OF RECORD 2 NOT REACHED.', 

* /IX,' NUMBER OF ENTRIES READ: ',14) 

4000 FORMAT( IX, ' A A* WARNING: OF RECORD 2 ENTRIES', 14, 

A ' NOT DIVISIBLE BY 7') 

END 


100 

200 

C 

225 
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SUBROUTINE READS YM < SYMM , LUNP > 

C 

CHARACTER A 80 RECORD 
REAL SYMM ( 3 ) 

LOGICAL LSYM1 ,LSYM2 ,LSYM3 
C 

REWIND LUNP I START SEARCH FROM TOP OF FILE 
C SET SYMMETRY DEFAULTS 
SYMM(l) =0.0 
SYMM( 2) a 0.0 
SYMM( 3 ) n 0.0 

C SEARCH PUNCH FILE FOR SYMMETRY DEFINITION RECORDS 
DO 100 I»l, 1000000 

READ( LUNP, 1000 ,END“150 ) RECORD 
IF (REC0RD(21:28) ,EQ. ' SYM1PLAN' ) THEN 
LSYM1 « .TRUE. 

DEC0DE( 2 , 2000 ,RECORD( 47 i 48) ) SYMM(l) 

TYPE A, 'PLANE 1 SYMMETRY : ' , SYMM(l) 

ELSE IF (RECORD< 21 i2Q) .EQ. ' SYM2PLAN ' ) THEN 
LSYM2 a .TRUE. 

DECODE ( 2, 2 000, RECORD f 47:40) > SYMM<2> 

TYPE A, 'PLANE 2 SYMMETRY:', SYMM(2) 

ELSE IF (RECORD! 21: 28) .EQ. ' SYM3PLAN' ) THEN 
LSYM3 = .TRUE. 

DECODE ( 2 , 2000 , RECORD ( 47:48) ) SYMM(3) 

TYPE A, 'PLANE 3 SYMMETRY:', SYMM<3) 

END IF 

100 CONTINUE 

C REPORT NON-DEFINED SYMMETRY 

150 IF (.NOT. LSYM1 ) TYPE 3000, '1' 

IF (.NOT. LSYM2) TYPE 3000, '2' 

IF (.NOT. LSYM3 ) TYPE 3000, '3' 

RETURN 

C 

1000 FORMAT (A80) 

2000 FORMAT (F2.0) 

3000 FORMAT( IX, 'PLANE ',A1,' SYMMETRY NOT SPECIFIED') 
END 




SUBROUTINE SORT !N, FREQ, GMASS, STIFF) 

C 

REAL FREQ(N) ,GMASS(N) ,STIFF(N) 

C 

C PERFORM ASCENDING SORT OF FREQ VECTOR 
C SLAVE SORT GMASS AND STIFF VECTORS 
C 

IF <N .EQ. 1) RETURN 
C 

DO 100 1=1, N-l 
DO 100 K°I+1 ,N 

IF ( FREQ ( K ) .LT. FREQ(I)> THEN 
TEMP = FREQ ( I ) 

FREQ ( I ) = FREQ(K) 

FREQ(K) « TEMP 
TEMP o GMASS ( I ) 

GMASS ( I) = GMASS (K) 

GMASS <K) = TEMP 
TEMP = STIFF ( I ) 

STIFF! I) = STIFF(K) 

STIFF(K) = TEMP 
END IF 

100 CONTINUE 

RETURN 
END 


SUBROUTINE WRITfiMTX ( N ,S?’REQ , GMASS , STIFF , SYMM , LUNM , LAMA > 
C 

CHARACTER MFILEMO , LAMA*8 

REAL FREQ(N) ,GMASS(N) ,STIFF(N> ,SYMM(3) 

C 

C OPEN MATRIX FILE 

TYPE 1000, 'ENTER OUTPUT MATRIX FILENAMES' 

ACCEPT 2000, MFILE 

OPEN (UNIT=LUNM,NAME=MFILE,TYPE= 'NEW' , 

* CARR I AGECONTROL = 'LIST' ) 

C WRITE MATRIX FILE HEADER 

IROWS = N 
JCOLS = 6 

WRITE( LUNM, 3000) LAMA, IROWS , JCOLS 
C WRITE VECTORS AS THREE COLUMNS OF MATRIX 
WRITE( LUNM, 4000) FREQ II COLUMN 1 
WRITE! LUNM, 4000) GMASS It COLUMN 2 
WRITE! LUNM, 4000) STIFF II COLUMN 3 
WRITE (LUNM, 40 00 ) ( SYMM! 1 ), 1=1 ,N) ll COLUMN 4 
WRITE! LUNM, 4000) (SYMM! 2 > ,1=1 ,N) I l COLUMNS 
WRITE! LUNM, 4000) ( SYMM! 3 ) , 1=1 ,N) I l COLUMN 6 
CLOSE! UNIT=LUNM) 

RETURN 

C 

1000 FORMAT! IX, A30, IX, $) 

2000 FORMAT (A40) 

3000 FORMAT! AO, 218) 

4000 FORMAT! 1PE16. 8) 

END 
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APPENDIX F. 


FORTRAN SOURCE CODE LISTING ( UNPACKDMI . FOR ) 


PROGRAM UNPACKDMI 

PARAMETER (MAXMTXSIZ «= lOOOOO) 

CHARACTER MATNAMA8, FLAGA1, CONTID* 3 
REAL ARRAY (MAXMTXSIZ) 

LOGICAL ALL, EOF 
DATA LUN / 10 / 

DO 500 NFILE°1 ,1000000 
GET NEXT FILENAME 
CALL GETFIL(LUN,*499) 

ALL « .FALSE. 

EOF = .FALSE. 

DO 400 NMATRX=1 , 1000000 

GET NEXT MATRIX NAME 
IF (.NOT. ALL) THEN 
TYPE 1000 

ACCEPT 2000, MATNAM 

IF (MATNAM .EQ. ' *ALL ') ALL = .TRUE. 

END IF 

SEARCH FILE FOR MATRIX AND GET DIMENSIONS 
CALL FNDMTX ( MATNAM , MROWS , NCOLS , CONTID , LUN , ALL , EOF ) 

IF (EOF) GOTO 399 

IF (MROWS*NCOLS .GT. MAXMTXSIZ) THEN 

TYPE A, 'MATRIX IS TOO LARGE', MATNAM, MROWS , NCOLS 
TYPE * , ' INCREASE PARAMETER MAXMTXSIZ IN UNPACKDMI' 
GOTO 399 
ENDIF 

UNPACK MATRIX INTO ARRAY 

CALL UNPACK (ARRAY , MROWS , NCOLS , CONTID , LUN , ALL , EOF ) 
REDUCE SQUARE DIAGONAL MATRIX TO VECTOR ? 

IF (MROWS .EQ. NCOLS) THEN 
TYPE 6000, MATNAM 
ACCEPT 4000, FLAG 
IF (FLAG .EQ. 'Y') 

A CALL REDUCE (MATNAM, ARRAY, MROWS, NCOLS) 

END IF 

WRITE UNPACKED MATRIX TO NEW FILE 
CALL WRTMTX ( MATNAM , ARRAY , MROWS , NCOLS > 

IS ANOTHER MATRIX IN THIS FILE DESIRED ? 

IF (ALL .AND. EOF) GOTO 450 
IF (.NOT. ALL) THh$ 

TYPE 3000 
ACCEPT 4000, FLAG 
IF (FLAG .NE. 'Y') GOTO 450 
END IF 
CONTINUE 

NO MORE MATRICES DESIRED FROM THIS FILE 
CLOSE (UNIT=LUN) 


C IS ANOTHER FILE DESIRED ? 

499 TYPE 5000 

ACCEPT 4000, FLAG 

IF (FLAG ,NE. ' Y ' ) GOTO 550 

500 CONTINUE 

550 STOP 'NO MORE FILES REQUESTED. ' 

C 

1000 FORMAT ( ' ENTER MATRIX NAME. ("*ALL" FOR ALL MATRICES): ',$) 

3000 FORMAT (AS ) 

3000 FORMAT < ' DO YOU WANT ANOTHER MATRIX IN THIS PILE ? (Y OR B)i ',$) 
4000 FORMAT(Al) 

5000 FORMAT ( ' DO YOU WANT ANOTHER FILE ? (Y OR N> t ',$> 

6000 FORMAT ( ' REDUCE DIAG MATRIX ',A8,' TO A VECTOR ? (7 OR N> : ',$> 
END 


SUBROUTINE FNDMTX ( MATNAM , MROWS , NCOLS , CONTID , LUN , ALL , EOF ) 

C 

CHARACTER MATNAM *8 , RECORD* 80 , FLD1T03*24, CONTID*3 
LOGICAL ALL, EOF 
C 

IF (.NOT. ALL) THEN 
REWIND LUN 
EOF = .FALSE. 

END IF 
C 

C SEARCH FOR DESIRED MATRIX AND GET DIMENSIONS 
FLD1T03 = 'DMI '//MATNAM//' O' 

DO 100 1=1,1000000 

READ < LUN, 1 000 ,END= 199 ) RECORD 
IF (ALL) FLD1T03 (9:16) = RECORD(9:16) 

IF ( RECORD ( 1:24) .EQ. FLD1T03 ) THEN 
IF (ALL) MATNAM = REC0RD(9:16) 

READ (RECORD, 2 000) MROWS, NCOLS, CONTID 
CONTID { 1:1) = 'A' 

GOTO 150 
END IF 

100 CONTINUE 

150 RETURN 

C 

C END OF FILE? NEXT MATRIX NOT FOUND 
199 EOF = .TRUE. 

IF (ALL) TYPE *, 'END OF FILE REACHED' 

IF (.NOT. ALL) TYPE *, MATNAM, 'MATRIX HEADER CARD NOT FOUND' 
RETURN 
C 

1000 FORMAT (ABO) 

2000 FORMAT(56X,2I8,A3 ) 

END 
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SUBROUTINE GETFIL{ LUN , * > 

C 

CHARACTER FILNAMMO 
C 

C GET NEXT FILENAME 
TYPE 1000 

ACCEPT 2000, FILNAM 
C OPEN FILE 

0 PEN < U N IT* LUN , NAME^FI LHAM , TYPE =' OLD ' , READONLY, ERR=199 ) 
RETURN 
C 

C NEW DMI FILE NOT FOUND 

199 TYPE A, FILNAM, 'FILE NOT FOUND' 

RETURN i 
C 

1000 FORMAT ( ' ENTER NAME OF NEXT DMI FILEs ',$> 

2000 FORMAT (A40) 

END 


SUBROUTINE REDUCE t MATNAM , MATRIX , MROWS , NCOLS ) 

C 

CHARACTER A 8 MATNAM 
REAL MATRIX (MROWS, NCOLS) 

C 

C CHECK FOR SQUARE MATRIX 
IF (MROWS .NE. NCOLS) THEN 

TYPE A, MATNAM, 'NOT A SQUARE MATRIX. NO REDUCTION DONE.' 
RETURN 
ENDIF 
C 

C REDUCE TO (IX NCOL) ROW VECTOR USING DIAGONAL VALUES 
DO 100 N=l, NCOLS 

MATRIX(N,1) = MATRIX ( N,N) 

100 CONTINUE 

MROWS = 1 
RETURN 
END 
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SUBROUTINE UNPACK ( MATR IX , MRO WS , NCOLS , CONTID , LUN , ALL , EOF ) 
CHARACTER* 7 2 RECORD 

CHARACTER *16 FLDID(2i5), CARDID*B , C0NTID*3 
EQUIVALENCE ( RECORD! 1»1), CARDID(ltl) ) 

EQUIVALENCE ( REC0RD(9:9), FLDID! 2 ) { 1 ; 1 ) ) 

CHARACTER*6 FMT<6) , REAL, INTGR, BLANK 
REAL MATR IX ( MRO WS, NCOLS) 

LOGICAL ALL, EOF 

DIMENSION VALUE<4), IVALUE(4) 

EQUIVALENCE < VALUE! 1 ) , IVAL!JE( 1 ) ) 

DATA FMT ( 1 ) / ' < 8X , ' / , FMT! 6 ) / ' 8X) ' / 

DATA REAL/ 'E16. B, ' / , INTGR/ ' 116 ,' / , BLANK/ ' A16 ,' / 


ZERO THE MATRIX 
DO 100 1=1 ,MR0WS 
DO 100 J=l, NCOLS 

MATRIX! I , J) =0.0 
CONTINUE 


PROCESS THE DMI CARDS FOR THIS MATRIX 
DO 400 K=l, 1000000 

READ! LUN, 1000, END=425) CARDID,FLDID 
IF (CARDID .EQ. 'DMI ') THEN 

THIS IS HEADER CARD OF NEXT MATRIX; THIS MATRIX DONE. 
IF (ALL) BACKSPACE LUN 
GOTO 450 

ELSE IF (CARDID .EO. 'DMI* ') THEN 

THIS IS A COLUMN IDENTIFIER CARD 
READ ( RECORD , 2 0 0 0 ) JCOL , IROW , VALUE ( 1 ) 

MATRIX! IROW,JCOL) = VALUE! 1) 

IROW = IROW + 1 

ELSE IF (CARDID! li 3) .EQ. CONTID) THEN 
THIS IS A ROW/VALUE CARD 
DETERMINE TYPE OF CONTENTS IN FIELDS 2;5 
DO 200 L=2 , 5 

IF (FLDID(L) ( 13:13) .EQ. 'EM THEN 
FMT(L) = REAL 

ELSE IF (FLDID(L) (16:16) ,E0. ' ') THEN 
FMT(L) = BLANK 

ELSE 

FMT(L) = INTGR 
ENDIF 

200 CONTINUE 

C READ CARD IN APPROPRIATE FORMAT 

READ! RECORD,FMT,ERR=201) VALUE 
C PROCESS FIELDS 2 : 5 

201 DO 300 L=2 , 5 

IF (FMT(L) .EQ. REAL) THEN 

MATRIX! IROW, JCOL) = VALUE(L-l) 

IROW = IROW + 1 

ELSE IF ( FMT ! L ) .EQ. INTGR) THEN 
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IROW = IVALUE(L-l) 

ELSE IF (FMT(L) ,EQ. BLANK) THEN 
CONTINUE 
ENDIF 

300 CONTINUE 

ELSE 

C THIS IS A NON-DMI CARD; THIS MATRIX DONE 

READ ( RECORD ,1000) CARDID , FLDID 
GOTO 450 
ENDIF 

400 CONTINUE 

C END OF FILE 

425 EOF « .TRUE. 

C END OF MATRIX 

450 RETURN 

C 

1000 FORMAT (AS, 4 A16) 

2000 FORMAT ( 24X,2I16,E16.8) 

END 


SUBROUTINE WRTMTX ( MATNAM , MATRIX , MROWS , NCOLS ) 

CHARACTER MATNAM'*' 8, NAME>*<12 
REAL MATRIX (MROWS, NCOLS) 

C 

NAME » MATNAM// ' .MTX' 

OPEN ( UN IT= 5 0 , NAME=NAME , TYPE= ' NEW ' ,CARRIAGECONTROL= ' LIST ' ) 
C WRITE HEADER RECORD TO ASCII FILE 
WRITE ( 50,1000) MATNAM, MROWS, NCOLS 
C WRITE MATRIX TO ASCII FILE ( 1 VALUE PER RECORD ) 

WRITE( 50 , 2000 ) MATRIX 
CLOSE ( UNIT=50 ) 

TYPE *, 'MATRIX FILE WRITTEN WITH MROWS, NCOLS:' 

TYPE A, NAME, MROWS, NCOLS 
RETURN 
C 

1000 FORMAT ( A8 , 2 I 8 ) 

2000 FORMAT ( 1PE16 . 8 ) 

END 
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APPENDIX G. FORTRAN SOURCE CODE LISTING (GRDPTLST.FOR) 


PROGRAM GRDPTLST 
C 

C THIS PROGRAM WILL READ THE TESET (DOF ID VS. GRID POINT ID) 

C MATRIX FILE OBTAINED FROM THE NASTRAN PUNCH FILE AND GENERATE 
C A GRID POINT ID LIST FILE, USING THE NON-ZERO ENTRIES OF THE 
C TESET MATRIX. 

C 

CHARACTER ’“'35 MTXFIL, GPFIL , MTXNAME*8 
C 

C GET TESET MATRIX FILENAME 

10 TYPE 1000, 'ENTER GRID POINT ID MATRIX FILENAME: ' 

ACCEPT 2000, MTXFIL 

0PEN( UNIT=1 ,NAME=MTXFIL ,TYPE* ' OLD ' , READONLY , ERR =10 ) 

READ (1,3000) MTXNAME , MROWS, NCOLS 

IF (NCOLS .NE. 1) STOP 'ERROR: SHOULD BE ONLY ONE COLUMN' 

C GET GRID POINT LIST FILENAME 

TYPE 1000, 'ENTER GRID POINT ID OUTPUT LIST FILENAME:' 

ACCEPT 2000, GPFIL 

OPEN(UNIT=2,NAME=GPFIL, TYPE* 'NEW' , CARR IAGECONTRO L= ' LIST' ) 

C GENERATE LIST FILE FROM NON-ZERO MATRIX ENTRIES 
NIDS = 0 

DO 100 MD0F=1, MROWS 

READ( 1,5000) GRDPT 
IF (GRDPT .NE. 0.0) THEN 
IGRDPT = I NT (GRDPT) 

COMP = GRDPT - FLOAT (IGRDPT) 

ICOMP <= NINT(10.0*COMP) 

IF (ICOMP .GT. 6) THEN 
IGRDPT « IGRDPT + 1 
ICOMP = 0 
END IF 

WRITE( 2 ,6000 ) IGRDPT, ICOMP, MDOF 
NIDS = NIDS + 1 
END IF 

100 CONTINUE 

CL0SE(UNIT=1 ) 

CLOSE) UNIT=2) 

TYPE NIDS, 'ENTRIES WRITTEN TO GRID POINT LIST' 

STOP 

C 

1000 FORMAT ( IX, A44, IX, $) 

2000 FORMAT (A3 5 ) 

3000 FORMAT) AS, 218) 

4000 FORMAT ( IX, A8 , ' IS A '.IQ,' ROW BY ',18,' COLUMN MATRIX') 

5000 F0RMAT(E16.S) 

6000 FORMATdB,'-' ,11,112) 

END 
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APPENDIX H. FORTRAN SOURCE CODE LISTING ( STATCORR . FOR ) 


PROGRAM STATCORR 
C 

PARAMETER { MAXMODE*: 200) 

PARAMETER (MAXDOF » 1000) 

COMMON /LIMITS/ MXMD,MXDF 
DATA MXMD / MAXMODE / , MXDF/MAXDOF/ 

C 

REAL FREQ AN ( MAXMODE ) , FREQEX ( MAXMODE ) 

REAL DI SPAN {MAXDOF, MAXMODE) , DISPEX( MAXDOF, MAXMODE) 

REAL MASS (MAXMODE) , STIFF (MAXMODE) , DAMP (MAXMODE) 

REAL CORREL (MAXMODE, MAXMODE ) , C ( MAXMODE , MAXMODE ) 

REAL S( MAXMODE, MAXMODE) 

REAL RMSA{ MAXMODE, MAXMODE) , RMSE( MAXMODE, MAXMODE) 

REAL AVECT< MAXDOF) ,EVECT( MAXDOF) ,DIFF( MAXDOF) ,DIFABS (MAXDOF) 
INTEGER IDDOFA< MAXDOF, MAXMODE) ,IDD0FE( MAXDOF, MAXMODE) 

INTEGER ICOMPA (MAXDOF, MAXMODE ) , IC0MPE( MAXDOF, MAXMODE) 

INTEGER NDOFA< MAXMODE ) , NDOFEf MAXMODE) 

INTEGER IDDOFI ( MAXDOF ) ,ICOMP I (MAXDOF) 

INTEGER NAFIT ( MAXMODE ) , MEFIT ( MAXMODE ) 

INTEGER AS YMM <3, MAXMODE) , ES YMM ( 3 , MAXMODE ) 

CHARACTER’ 1 ' 10 AHEADER ( 4 ) , EHEADER < 4 ) , HEADER ( 3 2 > 

CHARACTER/' 9 CDATE 
CHARACTER *1 YESNO 
LOGICAL PRTAN , PRTEX , LS YMM , LPRINT 
C 

COMMON /WHEN/ CDATE 
C 

DATA PRTAN, PRTEX, LSYMM, LPRINT/ 4*. FALSE. / 

DATA AHEADER/ 2 *' NASTRAN ' , 2* ' TEST'/ 

DATA EHEADER/ 2 *' TEST', 2*' NASTRAN'/ 

DATA HEADER/' OFGP"S', ' CORREL', 2*' ',2*' RMS', 

A ' MAX REL' , ' GRID', ' MODE', 'FREQUENCY', 

A ' MODE', 'FREQUENCY', ' COMPARED', ' COEFF' , 

A ' CA', ' S', '(NASTRAN)', ' (TEST)', 

A 'DIFFERENCE', ' POINT', 12*' '/ 

C 

C PRINT SECTION 1 HEADER 

CALL DATE (CDATE) 

PRINT 8000, CDATE 
PRINT 4000 

C CHECK FOR SEPARATE PRINT FILE 

TYPE 5000 

ACCEPT 7000, YESNO 

IF (YESNO .EQ. 'Y'> LPRINT = .TRUE. 

IF (LPRINT) PRINT 5500, YESNO 

C GET ANALYTICAL FREQ, MASS, STIFFNESS, MODE-SHAPE, AND DOF LISTS 
CALL GETAN(FREQAN, MASS, STIFF, ASYMM, DISPAN, NAN, 

* IDDOFA, ICOMPA, NDOFA, LPRINT) 

C GET EXPERIMENTAL FREQ, DAMPING, MODE-SHAPE, AND DOF LISTS 
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CALL OETEXP ( FREQEX , DAMP , ES YMM , DISPEX , HEX , IDDOFE , ICOMPE , 

A NDOFE,LPRINT) 

C DETERMINE SYMMETRY, MODE-SHAPE DUMP, AND THRESHOLD OPTIONS 
CALL OPTIONS ( LS YMM , PRTAN , PRTEX , RTHRESH , LPRINT ) 

C PRINT INPUT SUMMARY 

CALL INPSUM ( NAN , FREQAN , MAS S , STIFF , AS YMM , MEX , FREQEX , DAMP ,ES YMM ) 

C SORT DOF LISTS AND MODE SHAPE MATRICES (BASED ON DOF SORT) 

CALL DOFSORT < IDDOFA , ICOMPA , NDOFA , DISPAN , NAN ) 

CALL DOFSORT ( IDDOFE , ICOMPE , NDOFE , DISPEX , MEX ) 

COMPUTE CORRELATION COEFFICIENTS AND RMS VALUES FOR ALL 

POSSIBLE ANALYTICAL /EXPERIMENTAL PAIRS 

DO 100 N=1 ,NAN 
DO 100 Mol, MEX 

GET VECTORS TO BE COMPARED AND THEIR DOF INTERSECTION SET 
CALL GETVEC( DISPAN (1,N) ,NDOFA( N) ,AVECT) 

CALL GETVECC DISPEX ( 1 ,M) , NDOFE ( M ) ,EVECT> 

CALL INTERSECT < AVECT , EVECT , IDDOFA <1,N) , IDDOFE U,M) , 

A ICOMPA ( 1 , N ) , ICOMPE { 1 , M ) , 

A NDOFA ( N ) ,NDOFE(M) ,IDDOFI ,IC0MPI ,NDOFI , 

A ASYMM(1,N) ,ESYMM(1,M) ,LSYMM) 

COMPUTE CORK COEFFS AND RMS VALUES FOR THIS ra AIR 
CALL CORRMS ( AVECT , EVECT , NDOFI , CORREL (N,M) , 

A C(N,M) ,S(N,M> ,RMSA(N,M) ,RM3E(N,M) ) 

i CONTINUE 

PRINT CORRELATION COEFFICIENT TABLE 
CALL CORRTBL( CORREL, NAN, MEX) 

DETERMINE BEST MATCH FOR EACH EXPERIMENTAL AND ANALYTICAL MODE 
CALL MATCH (CORREL, NAN, MEX, NAFIT,MEFIT) 

FOR EACH MATCHED PAIR, GET X/RMS DIFFERENCES AND PRINT SUMMARY 
PRINT ANALYTICAL MODE SHAPE HEADER 
PRINT 8000, CDATE 
PRINT 1000, AHEADER, HEADER 


DO 200 N=1 ,NAN 
M = NAFIT(N) 

C GET VECTORS TO BE COMPARED AND THEIR DOF INTERSECTION SET 

CALL GETVEC ( DISPAN (1,N) , NDOFA ( N ) , AVECT) 

CALL GETVEC ( DISPEX ( 1 ,K > , NDOFE (M) , EVECT ) 

CALL INTERSECT ( AVECT ,EVECT , IDDOFA ( 1 , N ) , IDDOFE ( 1 , M ) , 

* ICOMPA ( 1,N ) ,ICOMPE(l,M) , 

A NDOFA { N ) ,ND0FE(M> , IDDOFI , ICOMPI , NDOFI , 

* ASYMMf 1 ,N) ,ESYMM( 1 ,M) ,LSYMM) 

C CALCULATE INDIVIDUAL DIFFERENCES 

CALL RMSDIFF( AVECT, EVECT, NDOFI ,RMSA(N,M) ,RMSE(N,M) , 

A CORREL <N,M),S(N,M), DIFF , DIFABS ) 

C PRINT SUMMARY FOR THIS PAIR, WITH MAX DIFF AND > THRESHOLD 

CALL PRINT(N,FREQAN(N) ,M,FREQEX(M) ,ND0FI , CORREL (N,M) , 

A C(N,M) ,S(N,M) ,RMSA(N,M) ,RMSE(N,M) , 

A DIFF, DIFABS , IDDOFI , ICOMPI , RTHRESH) 


200 CONTINUE 

C PRINT EXPERIMENTAL MODE SHAPE HEADER 

PRINT 8000, CDATE 
PRINT 3000, EHEADER, HEADER 
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DO 300 M»1,MEX 
N » MEFIT(M) 

C GET VECTORS TO BE COMPARED AND THEIR DOF INTERSECTION SET 

CALL GETVEC(DISPEX( 1 ,M) ,NDOFE(M) , EVECT) 

CALL GETVEC ( DISPAN ( 1 , N ) ,NDOFA(N) , AVECT) 

CALL INTERSECT ( AVECT , EVECT , IDDOFA < 1 ,N > , IDDOFE ( 1 , M ) , 

A ICOMPA(l,N),ICOMPEU,M), 

A NDOFA(N) ,NDOFE(M) ,IDDOFI,ICOMPI ,ND0FI , 

A ASYMMf 1 ,N) ,ESYMM( 1 ,M> ,LSYMM) 

C CALCULATE COMBINED RMS AND INDIVIDUAL DIFFERENCES 

CALL RMSDIFF< AVECT, EVECT, NDOFI ,RMSA(N,M) ,RMSE<N,M) , 

A CORREL(N,M) ,S(N,M) ,DIFF,DIFABS > 

C PRINT SUMMARY FOR THIS PAIR, WITH MAX DIFF > THRESHOLD 

CALL PRINT (M,FREQEX(M) ,N,FREQAN(N) , NDOFI ,CORRELTn,M> , 

A C(N,M> ,S(N,M> ,RMSA(N,M) ,RMSE(N,M) , 

A DIFF ,DIFABS , IDDOFI , ICOMPI ,RTHRESH ) 

300 CONTINUE 

C PRINT ANALYTICAL AND/OR EXPERIMENTAL MODE-SHAPE VECTORS 
CALL PRINTPHI < NAN,NDOFA , IDDOFA , I C OMP A , D I S PAN , PRTAN , 

A MEX , NDOFE , IDDOFE , ICOMPE ,DISPEX , PRTEX > 

STOP 'FORTRAN STOP — PROCESSING COMPLETED' 

F0RMAT( 1H0 , ' 4 . ANALYTICAL MODE SHAPES AND THEIR BEST', 

A ' EXPERIMENTAL MATCHES? ' ,/'0' ,12(A10,1X> ,2(/12(lX,AI0) > ) 

FORMAT ( 1H0 , ' 5. EXPERIMENTAL MODE SHAPES AND THEIR BEST' 

A ' ANALYTICAL MATCHES}', /'O', 12(A10, IX), 2 < / 12 ( IX, A10 ) ) ) 

FORMAT { 1H0 , ' 1. INTERACTIVE DIALOG} ' ) 

FORMAT ( / , ' IS A SEPARATE OUTPUT LISTING FILE TO BE PRINTED?', 

A '(Y OR K)i ' ,$) 

FORMAT ( / , ' IS A SEPARATE OUTPUT LISTING FILE TO BE PRINTED?', 

A ' (Y OR N) S ' ,A1) 

FORMAT (Al) 

FORMAT ( 1H1 , 'NASTRAN MODAL ANALYSIS - MODAL SURVEY STATISTICAL', 

A ' CORRELATION', /1X,A9) 

END 


C 

1000 

3000 

4000 

5000 

5500 

7000 

8000 
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SUBROUTINE CORRMS < AVECT , EVECT , NDOFI , CORREL , CA , S , RMS A , RMSti ) 

COMPUTE CORRELATION COEFFICIENTS AND RMS VALUES 

REAL AVECT { A > , EVECT ( A ) 

INITIALIZE RMS'S AND CORR COEFFS 
RMSA » 0.0 
RMSE =0.0 
CORREL =0.0 
CA = 0.0 
S = 0.0 

IF (NDOFI .EQ. 0) RETURN 

COMPUTE VARIANCES OF ANALYTICAL AND EXPERIMENTAL VECTORS 
NOTE! THE FINAL VALUES COMPUTED HERE FOR THE "VARIANCES" 

ARE NOT TRUE VARIANCES. TO OBTAIN VARIANCES, THE VALUES 
VARA, VARE, AND VARAE SHOULD BE EACH MULTIPLIED BY 1 /NDOFI . 
THESE VALUES ARE NOT REPORTED, BUT ARE USED ONLY TO 
CALCULATE CORREL, CA, S, RMSA, AND RMSE. THE 1/ NDOFI FACTOR 
IS NOT INCLUDED, SINCE IT APPEARS IN BOTH THE NUMERATOR AND 
DENOMINATOR OF CORREL AND CA? THE EXTRA ARITHMETIC OPERATIONS 
WOULD SIMPLY INTRODUCE UNNECESSARY ERROR IN THE RESULTS. 

THE 1 /NDOFI FACTOR IS INCLUDED IN THE CALCULATIONS OF 
S, RMSA, AND RMSE. 

VARA =0.0 
VARE =0.0 
VARAE = 0.0 
DO 100 N=l, NDOFI 

VARA « VARA + AVECT ( N ) * AVECT ( N ) 

VARE = VARE + EVECT ( N ) * EVECT ( N ) 

VARAE = VARAE + AVECT(N)*EVECT(N) 

i continue; 

COMPUTE CORRELATION COEFFICIENT FOR THIS PAIR OF VECTORS 
IF (VARAE ,NE. 0.0) CORREL = VARAE/ SQRT ( VARA* VARE ) 

COMPUTE CORREL COEFF (REFERENCED TO ANALYTICAL VECTOR) 

IF (VARA .NE. 0.0) CA = VARAE/VARA 

COMPUTE ROOT OF MEAN SQUARE DIFFERENCE 
S = SQRT( ABS(VARE-CAACAAVARA)/ NDOFI ) 

COMPUTE RMS VALUES 

NOTE: FOR THIS APPLICATION, RMS = STD. DEV. = SQRT( VARIANCE) . 

RMSA = SQRT (VARA /NDOFI ) 

RMSE = SQRT (VARE /NDOFI) 

RETURN 

END 
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SUBROUTINE CORRTBL ( CORREL , NAN , MEX ) 

C 

C PRINT AN NAN X MEX TABLE OF THE ANALYTICAL VS. EXPERIMENTAL 
C CORRELATION COEFFICIENTS. 

C 

COMMON /LIMITS/ MAXMODE, MAXDOF 

COMMON /WHEN/ CDATE 

HEAL CORREL ( MAXMODE , MAXMODE ) 

CHARACTER CDATE/* 9 
C 

DO 200 MFIRST B 1 , MEX, 16 
DO 200 NFIRST° 1 , NAN ,50 
MLAST » MINI MEX, MFIRST+15 ) 

NLAST « MINI NAN, NFIRST+49 ) 

C PRINT HEADER FOR ONE PAGE OF THE TABLE 

PRINT 1000, CDATE 
PRINT 2000 

PRINT 3000, !M,M«MFIRST, MLAST) 

DO 100 N C NFIRST, NLAST 

PRINT 4000, N, ( CORREL ( N,M ) ,M=MFIRST, MLAST) 

100 END DO 

200 END DO 
C 

RETURN 

C 

1000 FORMAT! 1H1 , ' NASTRAN MODAL ANALYSIS - MODAL SURVEY STATISTICAL 
A 'CORRELATION', /1X,A9, 

A /1H0 , ' 3 . CORRELATION COEFFICIENTS FOR ANALYTICAL VS. ', 

A 'EXPERIMENTAL COMPARISONS! ' ) 

2000 FORMAT! ' OANALYTICAL ' ,24X, 'EXPERIMENTAL MODES') 

3000 FORMAT! 3X, 'MODES' ,3X,I5,15I7) 

4000 FORMAT! IX, 16, 4X, 16! IX, F6. 3) ) 

END 
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SUBROUTINE: D0F30RT( IDDOP , ICOMP , NDOF , DISP , NMODES ) 

PUT DOF LIST FOR EACH MODE INTO ASCENDING SORT, 

AND SLAVE-SORT THE CORRESPONDING MODE-SHAPE VECTORS. 

COMMON /LIMITS/ MAXMODE , MAXDOF 
REAL DISP (MAXDOF, MAXMODE) 

INTEGER IDDOF ( MAXDOF , MAXMODE) , I COMP (MAXDOF, MAXMODE) ,NDOF( MAXMODE) 
C 

DO 200 J»l, NMODES 

IF (NDOF(J) .EQ, 1) GOTO 200 
DO 100 I«1,ND0F(J>-1 
DO 100 KnI+l,NDOF(J> 

IF <IDD0F(K,J) ,LT. IDDOF(I,J)) THEN 
CALL SWAPC IDDOF( I , J) , IDDOF(K,J) ) 

CALL SWAP( ICOMPd , J) , IC0MP(K,J) ) 

CALL SWAP( DISP( I , J) , DISP(K,J) > 

ELSE IF ( ( IDDOF(K, J) ,EQ. IDDOF(I,J>) .AND. 

^ ( ICOMP(K, J) .LT. ICOMP(I, J)) > THEN 

CALL SWAP( IDDOFt I , J) , IDDOF(K,J) ) 

CALL SWAP( ICOMPd, J), ICOMP(K,J) ) 

CALL SWAP( DISP( I,J) , DISP(K,J) ) 

END IF 

100 CONTINUE 

200 CONTINUE 

RETURN 
END 
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SUBROUTINE GETAN ( FREOS , MASS , STIFF , ASYMM , SHAPES ,NM0DE3 f 
A IDDOF , ICOMPA , NDOFS , LPRINT ) 

C 

C GET FREQUENCY, MASS, STIFFNESS AND SYMMETRY VECTORS; 

C MODE-SHAPE MATRIX; AND GRID POINT LIST. 

C 

COMMON /LIMITS/ MAXMODE, MAXDOF 
C 

REAL FREQS { MAXMODE) ,MASS( MAXMODE) , STIFF! MAXMODE) 

REAL SHAPES (MAXDOF, MAXMODE) 

INTEGER ASYMM (3, MAXMODE) 

INTEGER IDDOF (MAXDOF, MAXMODE) , ICOMPA (MAXDOF, MAXMODE) 

INTEGER NDOFS (MAXMODE) 

CHARACTER A3 5 FRQFIL ,SHPFIL ,DOFFIL 
CHARACTER A 8 FRQMTX,SHPMTX 
LOGICAL LPRINT 
C 

C GET FILE NAMES 

10 TYPE 1000 

ACCEPT 2000, FRQFIL 

OPEN ( UNIT«1 , NAME=FRQFIL ,TYPE" ' OLD ' .READONLY, ERR°10 ) 

IF (LPRINT 1 ) PRINT 1500, FRQFIL 
20 TYPE 3000 

ACCEPT 2000, SHPFIL 

OPEN ( UNIT- a 2 ,NAME=SHPFIL ,TYPE e ' OLD' ,READ0NLY,ERR"20 ) 

IF (LPRINT) PRINT 3500, SHPFIL 
30 TYPE 4000 

ACCEPT 2000, DOFFIL 

OPEN ( UNIT» 3 , NAME=DOFFIL ,TYPE" ' OLD ' , READONLY , ERR *3 0 ) 

IF (LPRINT) PRINT 4500, DOFFIL 
C READ FREQUENCY FILE 

READ( 1,7000) FRQMTX,NFREQS ,NCOLS 
IF (NCOLS .NE. 6) THEN 
40 TYPE 5000, FRQMTX, NCOLS 

IF (LPRINT) PRINT 5000, FRQMTX, NCOLS 
CLOSE ( UNIT" 1) 

CLOSE (UNIT" 2) 

CL0SE(UNIT=3) 

STOP 
END IF 

IF (NFREQS ,GT. MAXMODE) PRINT 6000, NFREQS , MAXMODE 
NMODES " MIN (NFREQS, MAXMODE) 

DO 100 I«1 ,6 
DO 100 J=l, NFREQS 

IF (J .GT. KAXMODE) THEN 
READ( 1,8000) DUMMY 

ELSE IF (I .EQ. 1) THEN I GET FREQUENCY LIST 
READ! 1,8000) FREQS(J) 

ELSE IF (I .EQ. 2) THEN I GET MASS LIST 
READ( 1 ,8000 ) MASS(J) 

ELSE IF (I .EO. 3) THEN I GET STIFFNESS LIST 
READ( 1 ,8000 ) STIFF(J) 
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ELSE IF (I .EQ. 4) THEN I GET PLANE 1 SYMMETRY LIST 
READ( I / 8000 > SYMM 
ASYMM(1,J) “ NINT(SYMM) 

ELSE IF (I .EQ, 5) THEN I GET I’MtiJS 2 SYMMETRY LIST 
READf 1,8000) SYMM 
ASYMM(2,J> « NINT(SYMM) 

ELSE IF (I .EQ. 6) THEN ! GET PLANE 3 SYMMETRY LIST 
READf 1,8000) SYMM 
ASYMM(3,J) n NINT(SYMM) 

END IF 

100 CONTINUE 

C READ MODE-SHAPE FILE 

READ( 2 .7000 ) SHPMTX, NDISPS , NSHAPES 
IF (NSHAPES .NE. NFREQS ) THEN 

NMODES - MIN( NFREQS, NSHAPES, NMODES ) 

TYP F. 9000, NFREQS, NSHAPES, NMODES 

IF (LPRINT) PRINT 9000, NFREQS , NSHAPES, NMODES 

END IF 

IF (NDISPS , GT, MAXDOF) PRINT 10000, NDISPS ,MAXDOF 
DO 200 J»l, NMODES 
DO 200 lol, NDISPS 

IF (I .LE. MAXDOF) READ<2,8000> SHAPES ( I, J) 

IF (I .GT. MAXDOF) READf 2, 8000) DUMMY 
200 CONTINUE 

C READ DOF ID FILE 
NDOFSfl) ■ 0 
DO 300 Ini, 1000000 

IF (I .LE. MAXDOF) 

A READ( 3 , 11000, ENDo 350 ) IDDOF(I,l) ,ICOMPA(I,l) 

IF (I .GT. MAXDOF) READf 3 ,11000 ,END»350 ) IDUMMY, IDUMMC 
NDOFS(l) n NDOFS(l) + 1 
300 CONTINUE 

350 IF (NDOFS(l) .NE. NDISPS) THEN 

TYPE 12000, NDISPS, NDOFS(l) 

IF (LPRINT) PRINT 12000, NDISPS, NDOFSf 1 ) 

CLOSEf UNIT=1) 

CL0SE(UNIT=>2 ) 

CLOSEf UNIT=3 ) 

STOP 

ENDIF 

NDOFS(l) = MIN (NDOFSfl), MAXDOF) 

IF (NMODES .EQ. 1) RETURN 
DO 500 J=2, NMODES 

NDOFSf J) s NDOFS(l) 

DO 400 I=l,NDOFS( J) 

IDDOFf I ,J ) = IDDOFf 1,1) 

ICOMPAf I, J > = ICOMPA( 1,1 ) 

400 CONTINUE 

500 CONTINUE 

CLOSE( UNITED 
CLOSE(UNTT=2 ) 

CLOSEf UNIT 3 3 ) 
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RETURN 

C 

1000 FORMAT* ENTER ANALYTICAL LAMA MATRIX FILENAME; ',$) 

1500 FORMAT*/,' ENTER ANALYTICAL LAMA MATRIX FILENAME; ' ,A35) 

2000 FORMAT* A3 5) 

3000 FORMAT*' ENTER ANALYTICAL MODE-SHAPE MATRIX FILENAME; ',$> 

3500 FORMAT*' ENTER ANALYTICAL MODE-SHAPE MATRIX FILENAME; ',A35) 

4000 FORMAT*' ENTER ANALYTICAL GRID POINT LIST FILENAME; ',$) 

4500 FORMAT*' ENTER ANALYTICAL GRID POINT LIST FILENAME; ' ,A35) 

5000 FORMAT*' FREQUENCY MATRIX ',A8,' HAS ',16,' COLUMNS.' 

A ' SHOULD BE 6 COLUMNS . ' ) 

6000 FORMAT*' AAAA WARNING; FREQUENCY VECTOR HAS ',13,' ENTRIES.', 

A /15X, 'ONLY THE FIRST', 15,' WILL BE USED.') 

7000 FORMAT* AB, 218) 

8000 FORMAT* E16. 8) 

9000 FORMAT*' AAAA WARNING; UNEQUAL NUMBER OF FREQUENCIES AND ' 

A 'MODE SHAPES. ', 

A /15X, 'NUMBER OF FREQUENCIES ;' ,16 , 

A /15X, 'NUMBER OF MODE SHAPES; ',16, 

A / 15X, ' ONLY THE FIRST', 15,' WILL BE USED.') 

10000 FORMAT*' AAAA WARNING; MODE-SHAPE VECTORS HAVE ',13,' ENTRIES.', 
A / 15X, ' ONLY THE FIRST', 16,' WILL BE USED.') 

11000 FORMAT* 18 , IX, II ) 

12000 FORMAT* ' AAAA ERROR; UNEQUAL NUMBER OF MODE-SHAPE POINTS AND ', 

A 'D.O.F. ID"S. ',16,' ',16) 

END 
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SUBROUTINE GETEXP(FREQE< .DAMP ,ESYMM , DISPEX , MEX , 

A IDDOFE, IC0MP,ND0FE, LPRINT) 

C 

C GET EXPERIMENTAL FREQUENCY DAMPING LISTS, MODE-SHAPE MATRIX, 

C AND D.O.F. LIST. 

C 

COMMON /LIMITS/ MAXMODE , MAXDOF 
C 

REAL FREQEX( MAXMODE) , DAMP ( MAXMODE ) , DISPEX (MAXDOF, MAXMODE) 

REAL DUMMY( 5) 

INTEGER ESYMMC 3, MAXMODE) 

INTEGER I DDOFE( MAXDOF, MAXMODE ) ,ICOMP< MAXDOF, MAXMODE) 

INTEGER NDOFE( MAXMODE) 

CHARACTER A 3 5 MSFILE 
LOGICAL LPRINT 
C 

C PROCESS MODE-SHAPE FILES 
TYPE *, ' ' I SKIP A LINE 
IF (LPRINT) PRINT A, ' ' 

MEX = 0 

DO 200 J=l, MAXMODE 

C GET NEXT EXPERIMENTAL MODE-SHAPE FILE 

10 TYPE 1000 

ACCEPT 2000, MSFILE 

IF (MSFILE ,EQ. 'NONE') GOTO 250 

OPEN ( UNIT* 1 0 , NAME= MSFILE , TYPE= ' OLD ' , READONLY , ERR = 1 0 ) 

IF (LPRINT) PRINT 1500, MSFILE 

C GET FREQUENCY, DAMPING AND SYMMETRY FOR THIS MODE-SHAPE 

READ110,*) FREQEX(J) 

READUO,*) DAMP ( J ) 

READ( 10 , A } (ESYMM(K,J) ,K=1,3) 

C GET MODE-SHAPE DISPLACEMENTS AND D.O.F. 'S 

DO 100 1=1, MAXDOF 

READ(10,*,END=150> 

* IDDOFEt I , J) , ICOMP ( I , J > , DISPEX ( I , J) 

NDOFE(J) * I 

100 CONTINUE 

PRINT 5000, MAXDOF, MAXDOF 
150 MEX = J 

CLOSE( UNIT=10 ) 

200 CONTINUE 

TYPE 6000, MAXMODE 
IF (LPRINT) PRINT 6000, MAXMODE 

250 IF (LPRINT .AND. (MSFILE ,EQ. 'NONE')) PRINT 1500, MSFILE 
IF (MEX .EQ. 0) THEN 

TYPE *,' AAA A ERROR: NO EXPERIMENTAL MODE-SHAPE REQUESTS' 

IF (LPRINT) PRINT A, 

A ' AAAA ERROR: NO EXPERIMENTAL MODE-SHAPE REQUESTS' 

STOP 
END IF 
RETURN 


1000 FORMAT (' ENTER NEXT EXPERIMENTAL MODE-SHAPE FILENAME. ', 

A ' ( "NONE" IF NO MORE) j ' ,$> 

1500 FORMAT ( ' ENTER NEXT EXPERIMENTAL MODE-SHAPE FILENAME. 

A '("NONE" IF NO MORE) » ',A35> 

2000 FORMAT (A35) 

5000 FORMAT ( ' AAAA WARNING! MODE-SHAPE VECTOR HAS' ,16,' OR MORE', 

A ' ENTRIES.' ,/13X,' ONLY THE FIRST', 16,' WILL BE USED.') 

6000 FORMAT ( ' AAAA NO MORE EXPERIMENTAL MODE-SHAPE FILES PERMITTED. 
A ' MAXIMUM IS ' ,14 > 

END 
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SUBROUTINE GETVEC(DISP,NDOF,VECT) i 

C ( 

C COPY VECTOR OF LENGTH NDOF FROM PERMANENT VECTOR TO 
C WORKING VECTOR. 

C 

REAL DISP(A) ,VECT( A) j 

C t 

IF (NDOF .EQ. 0) RETURN 

DO 100 1=1, NDOF 

VECT(I) = DISP(I) 

100 CONTINUE 

RETURN 

END 
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SUBROUTINE INPSUM( NAN , FREQAN , HASS , STIFF , ASYMM , 

A MEX , FREQEX , DAMP , ESYMM ) 

C 

C PRINT A SUMMARY OF THE ANALYTICAL AND EXPERIMENTAL FREQUENCIES, 
C MASS, STIFFNESS, DAMPING, AND SYMMETRY VECTORS. 

C 

REAL FREQAN <M ,MASS<*) , STIFF! A) 

REAL FREQEX(A) , DAMP (A) 

INTEGER AS YMM ( 3 , A ) , E5 YMM ( 3 , A ) 

LOGICAL NOSYMA<3> ,YESYMA(3) ,N0SYME(3> ,YESYME<3) .FIRST 
CHARACTER CDATE*9 
COMMON /WHEN/ CDATE 

DATA NOSYMA, YESYMA,NOSYME,YEKYME/12A. FALSE. /, FIRST/ .TRUE. / 

C 

C PRINT HEADER 

PRINT 8000, CDATE 
PRINT 1000 

C PRINT ANALYTICAL SUMMARY 

C PRINT ANALYTICAL HEADER 

PRINT 2000 

C PRINT SUMMARY DATA 

PRINT 3000, ( I , FREQAN ( I) ,MASS( I) , STIFF ( I > , 

A (ASYMM( J,I) ,J=1,3) , 1=1, NAN ) 

C CHECK FOR INCONSISTENT ANALYTICAL "UNDEFINED" SYMMETRY 

DO 200 IAXS=1 ,3 

DO 100 MODEa l r NAN 

IF ( ASYMM< IAXS ,MODE) ,EQ. 0 ) NOS YMA( IAXS )= .TRUE. 

IF ( AS YMM (IAXS, MODE ) .NE.O) YES YMA( IAXS) * .TRUE. 

100 CONTINUE 

IF ( NOSYMAt IAXS) .AND. YESYMA( IAXS) > PRINT 6000, IAXS 
200 CONTINUE 

C PRINT EXPERIMENTAL SUMMARY 
C PRINT EXPERIMENTAL HEADER 

PRINT 4000 

C PRINT SUMMARY DATA 

PRINT 5000, ( I ,FREQEX( I ) ,DAMP< I ) , 

A ( ESYMM (J,I) ,J=1,3) , 1=1, MEX) 

C CHECK FOR INCONSISTENT EXPERIMENTAL "UNDEFINED" SYMMETRY 

DO 400 IAXS=1 , 3 

DO 300 M0DE=1 ,MEX 

IF ( ES YMM (IAXS, MODE ) .EQ.O) NOS YME( IAXS )=. TRUE. 

IF ( ESYMM (IAXS, MODE) .NE.O) YES YME ( IAXS )=. TRUE . 

300 CONTINUE 

IF ( NOS YME ( IAXS) .AND. YESYME( IAXS) ) PRINT 6000, IAXS 
400 CONTINUE 

C CHECK FOR INCONSISTENT ANA VS. EXP "UNDEFINED" SYMMETRY 
DO 500 IAXS=1 ,3 

IF ( (NOSYMA (IAXS) .AND. YES YME ( IAXS ) ) .OR. 

A ( YESYMA (IAXS). AND. N0SYME( IAXS > ) ) THEN 

IF (FIRST) PRINT 7000 
FIRST = .FALSE. 

PRINT 6000, IAXS 


500 

c 

1000 

2000 


3000 

4000 


5000 

6000 

7000 

0000 


END IF 

CONTINUE 

RETURN 

FORMATQHO, '2. SUMMARY OF FREQUENCY, MASS, STIFFNESS , ' , 

* ' DAMPING, AND SYMMETRY:') 

FORMAT < / / , 5X , ' ANALYTICAL MODES : ' 

* /1H0,4X, 'MODE FREQUENCY MASS 

* ' STIFFNESS ' , ' SYMMETRY 

* / r 5X, ' ' ,4(1X, ')) 


FORMAT ( 
FORMAT ( 


A 

A 


FORMAT < 
FORMAT ( 

tc 

FORMAT ( 
FORMAT ( 
V 

END 


<NAN></5X,I3,1X,3UPE16.8) ,3X,3I3> ,/> 

II, 5X, ' EXPERIMENTAL MODES : ' , 

/ 1H0 , 4X , ' MODE FREQUENCY DAMPING SYMMETRY ' , 

/ # 5X , ' ' > 3 ( IX , ' ' ) ) 

<MEX>(/5X7l3, , lX,2aPE11.3) , IX, 313) ,/) 

10X, ' AAAWARNING: INCONSISTENT UNDEFINED SYMMETRY', 

' FOR PLANE' ,12) 

// ,5X, 'ANALYTICAL VS. EXPERIMENTAL SYMMETRY:') 

1H1 , 'NASTRAN MODAL ANALYSIS - MODAL SURVEY STATISTICAL 
'CORRELATION', /IX, A9) 


no 


a 


ORIGINAL F'A^t I2f 
OF POOR QUALTsf 


SUBROUTINE INTERSECT < AVECT , EVECT , IDDOFA , IDDOFE , ICOMPA , ICOMPE , 

A NDOFA , NDOFE , IDDOFI , ICOMPI , NDOFI , 

A AS YMM ,ESYMM , LSYMM ) 

C 

C COPY INTERSECTION SET OF VECTORS IDDOFA AND IDDOFE TO VECTOR 
C IDDOFI. NDOFI IS THE LENGTH OF THE INTERSECTION SET. 

C REDUCE AVECT AND EVECT VECTORS TO CONTAIN ONLY THE CORRESPONDING 

C INTERSECTION SET OF TERMS, I.A.W. THE MASTER SET IDDOFI, 

C IF THE INTERSECTION SET IS NULL (NO EQUIVALENT D.O.F ID'S), 

C THE VALUE RETURNED FOR NDOFI IS 0 . 

C IF SYMMETRY IS TO BE CONSIDERED, AND ANALYTICAL AND EXPERIMENTAL 

C SYMMETRY ARE NOT THE SAME, THE VALUE RETURNED FOR NDOFI IS 0 . 

C 

REAL AVECT ( A ) , EVECT ( A ) 

INTEGER IDDOFA ( A ) , IDDOFE ( A ) , IDDOFI ( A ) 

INTEGER ICOMPA ( A ) , ICOMPE( A ) , ICOMPI ( A ) 

INTEGER ASYMM( 3 ) ,ESYMM< 3 ) 

LOGICAL LSYMM 

IA = 1 
IE = 1 
NDOFI = 0 
IF (LSYMM) THEN 

SYMMETRY IS TO BE CONSIDERED t 

IF A AND E ARE UNLIKE, LEAVE INTERSECTION SET AS NULL 
IF (ASYMM(l) .NE. ESYMMfl)) RETURN 
IF (ASYMM< 2) .NE. ESYMM( 2 ) ) RETURN 
IF (AS YMM ( 3 ) .NE. ESYMM(3)> RETURN 
END IF 

DO 100 1=1 ,NDOFA+NDOFE 

IF ( (IA.GT. NDOFA) .OR. ( IE. GT. NDOFE) ) GOTO 150 
IF ( ( IDDOFA ( IA) .EQ. IDDOFE ( IE) ) .AND. 

A ( ICOMPA (IA) .EQ. ICOMPE (IE)) ) THEN 

NDOFI = NDOFI + 1 
IDDOFI (NDOFI) = IDDOFA (IA) 

ICOMPI (NDOFI) = ICOMPA (IA) 

AVECT (NDOFI) = AVECT (IA) 

EVECT (NDOFI) = EVECT (IE) 

IA » IA + 1 
IE = IE +1 

ELSE IF ( IDDOFA (IA) .LT. IDDOFE (IE)) THEN 
IA = IA + 1 

ELSE IF (IDDOFE (IB) .LT. IDDOFA ( IA) ) THEN 
IE ■ IE + 1 
END IF 
100 CONTINUE 

150 RETURN 
END 
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SUBROUTINE HATCH ( CORREL , NAN , MEX , NAFIT , MEFIT > 

C 

COMMON /LIMITS/ MAXMODE , MAXDOF 
C 

DIMENSION CORREL (MAXMODE, MAXMODE) , NAF IT C MAXMODE) ,MEFIT( MAXMODE) 
C 

C DETERMINE BEST EXPERIMENTAL MODE-SHAPE MATCH (HIGHEST 
C CORRELATION COEFFICIENT) FOR EACH ANALYTICAL MODE-SHAPE. 

C 

DO 200 N"1,NAN 
BEST ■ -1.0 
DO 100 M=1 ,MEX 

IF (ABS ( CORREL < N,M ) > ,GT. BEST) THEN 
BEST o ABS ( CORREL ( N,M) > 

NAFIT(N) = M 
END IF 

100 CONTINUE 

200 CONTINUE 

C 

C DETERMINE BEST ANALYTICAL MODE-SHAPE MATCH (HIGHEST 
C CORRELATION COEFFICIENT) FOR EACH EXPERIMENTAL MODE-SHAPE. 

C 

DO 400 M«1,MEX 
BEST * -1.0 
DO 300 N=1 ,NAN 

IF ( ABS ( CORREL ( N,M ) ) ,GT. BEST) THEN 
BEST = ABS ( CORREL ( N,M ) ) 

MEFIT(M) » N 
END IF 

300 CONTINUE 

400 CONTINUE 

RETURN 
END 
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SU BROUTINE OPTIONS (LSYMM, PRTAN , PRTEX , RTHRESH , LPR INT ) 

C 

C DETERMINE SYMMETRY, MODE-SHAPE DUMP, AND THRESHOLD OPTIONS 
LOGICAL LSYMM, PRTAN, PRTEX, LPRINT 
CHARACTER* 1 YESNO 

C DETERMINE IF SYMMETRY IS TO BE CONSIDERED 
TYPE 1000 

ACCEPT 2000, YESNO 

IF (YESNO .EQ. 'Y') LSYMM = .TRUE. 

IF (LPRINT) PRINT 1500, YESNO 

C DETERMINE IF ANALYTICAL AND/OR EXPERIMENTAL MODE-SHAPE 
C VECTORS ARE TO BE PRINTED 

TYPE *, ' ' I SKIP A LINE 
IF (LPRINT) PRINT *, ' ' 

TYPE 3000, 'ANALYTICAL' 

ACCEPT 2000, YESNO 

IF (YESNO .EQ. ' Y' ) PRTAN = .TRUE. 

IF (LPRINT) PRINT 3500, 'ANALYTICAL', YESNO 
TYPE 3000, 'EXPERIMENTAL' 

ACCEPT 2000, YESNO 

IF (YESNO .EQ. 'Y') PRTEX » .TRUE. 

IF (LPRINT) PRINT 3500, 'EXPERIMENTAL', YESNO 
C GET RELATIVE DEVIATION THRESHOLD. DEFAULT IS 5.0 % 

RTHRESH = .05 
TYPE 4000, RTHRESH 
ACCEPT *, RTHRESH 
TYPE 5000, RTHRESH 
IF (LPRINT) THEN 

PRINT 4500, .05, RTHRESH 
PRINT 5000, RTHRESH 
END IF 
RETURN 
C 

1000 FORMAT ( / , ' IS ANALYTICAL VS. EXPERIMENTAL SYMMETRY TO BE ' , 

* 'CONSIDERED? (Y OR N) : ',$) 

1500 FORMAT ( / , ' IS ANALYTICAL VS. EXPERIMENTAL SYMMETRY TO BE ', 

A 'CONSIDERED? (Y OR N>: ' ,A1) 

2000 FORMAT (Al) 

3000 FORMAT { ' PRINT ',A12,' MODE-SHAPE VECTORS? (Y OR N> s ',$> 
3500 FORMAT ( ' PRINT ',A12,' MODE-SHAPE VECTORS? (Y OR N) : ',A1> 
4000 FORMAT ( / / ' RELATIVE DEVIATIONS GREATER THAN A THRESHOLD', 

A ' VALUE WILL BE PRINTED. ' , 

A I' THE DEFAULT THRESHOLD IS ',F6.3, 

A /' ENTER DESIRED THRESHOLD. “ , " FOR DEFAULT: ',$> 

4500 FORMAT ( / / ' RELATIVE DEVIATIONS GREATER THAN A THRESHOLD', 

A ' VALUE WILL BE PRINTED. ' , 

A /' THE DEFAULT THRESHOLD IS ',F6.3, 

A /' ENTER DESIRED THRESHOLD. FOR DEFAULT: ',F6.3> 

5000 FORMAT ( ' OREL ATI VE DEVIATIONS GREATER THAN ',2PF7.2, 

A ' % WILL BE PRINTED. ' > 

END 
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SUBROUTINE PRINT(N1 ,FREQ1 ,N2,FREQ2 ,NDOF ,CORREL ,C ,S ,RMSA ,RMSE, 

* DIFF , DIFABS , IDDOF , ICOMP , RTHRE3H ) 

C 

C PRINT CORRELATION INFORMATION FOR MODES N1 ,N2 COMPARISON 
C 

DIMENSION DIFF(A) „DIFABS(A) , IDDOF ( * ) , ICOMP ( * ) 

DIMENSION ID( 30000 ) 

C 

C DETERMINE MAXIMUM DIFFERENCE AND DIFFS > THRESHOLD 
DIFMAX =0.0 
IDMAX a 0 
NDIFF a 0 

IF (NDOF .GT. 0) THEN 
DIFMAX a ABS(DIFFU)) 

IDMAX = IDDOF(l) 

DO 100 N=1 ,NDOF 

IF ( AB3 (DIFF(N) ) .GT. DIFMAX) THEN 
DIFMAX = ABS (DIFF(N) ) 

IDMAX a IDDOF (N) 

END IF 

IF (ABS(DIFF(N) > .GT. RTHRESH) THEN 
NDIFF « NDIFF + 1 
ID (NDIFF) a N 
END IF 

100 CONTINUE 

END IF 
C 

C PRINT SUMMARY LINE FOR THIS MATCHED PAIR 
PRINT 1000, N1,FREQ1,N2,FREP2,ND0F,C0RREL,C,S, 

A RMSA,RMSE, DIFMAX, IDMAX 

IF (NDOF .EQ. 0) RETURN 
C 

C PRINT ALL DIFFS > THRESHOLD , DEFAULT IS 5% 

IF (NDIFF .GT. 0) THEN 

PRINT 2000, RTHRESH, ( IDDOF ( ID (N)), ICOMPf ID(N) ) , 

A DIFF(ID(N) ), N=l, NDIFF) 

PRINT 3000, (IDDOFdD(N) ) , ICOMP(ID(N) ) , DIFABS (ID(N) ) , 

* N*l, NDIFF) 

END IF 

RETURN 

FORMAT ( '0' ,I6,2(F15.6,I7> ,4X,F9 . 3 , 2X,5( 1PE11, 3 ) ,19) 

FORMAT( 2 IX, 'RELATIVE DEVIATIONS (XA/RMSA-XE/RMSE) > ', 

A 2PF7 . 2 , ' % j (GRID ID/DEVIATION)', 

* 20(/21X,5(I9, ,11,'/' , 1PE10 . 3 ) ) ) 

FORMAT(21X,' SCALED DIFFERENCES ( <XA-XE)/S ) > THRESHOLD s', 

A ' (GRID ID/DIFFERENCE) ' 

A 20(/21X,5d9,'-' ,11,'/' ,1PE10.3)> ) 

END 


C 

1000 

2000 


3000 


SUBROUTINE PRINTPHI ( NAN , NDOFA , IDDOFA , ICOMPA, DISPAN, PRTAN , 

A MEX , NDOFE , IDDOFE , ICOMPE , DISPEX , PRTEX ) 

C 

COMMON /LIMITS/ MAXMODE, MAXDOF 
COMMON /WHEN/ CDATE 
C 

REAL DI SPAN (MAXDOF, MAXMODE) ,DISPEX( MAXDOF, MAXMODE) 

INTEGER IDDOFA ( MAXDOF, MAXMODE ) , IDDOFE ( MAXDOF , MAXMODE ) 

INTEGER ICOMPA (MAXDOF, MAXMODE ) , ICOMPE (MAXDOF, MAXMODE) 

INTEGER NDOFA ( MAXMODE ) , NDOFE ( MAXMODE ) 

LOGICAL PRTAN, PRTEX 
CHARACTER CDATE* 9 
C 

IF (PRTAN) THEN I PRINT ANALYTICAL MODE-SHAPE VECTORS 
C PRINT HEADER 

PRINT 5000, CDATE 
PRINT 1000 

C PRINT ALL VECTORS 

DO 100 1=1, NAN 

PRINT 2000, I , ( IDDOFA < J, I ) , ICOMPA ( J , I ) , DISPAN (J ,1 ) , 

* J=1 ,ND0FA( I ) ) 

100 CONTINUE 

END IF 
C 

IF (PRTEX) THEN I PRINT EXPERIMENTAL MODE-SHAPE VECTORS 

C PRINT HEADER 

PRINT 5000, CDATE 
PRINT 3000 

C PRINT ALL VECTORS 

DO 200 1=1, MEX 

PRINT 4000, I , ( IDDOFE (J, I > , ICOMPE ( J, I ) , DISPEX (J ,1 > , 

* J=1 ,NDOFE( I ) ) 

200 CONTINUE 

END IF 
RETURN 
C 

1000 FORMAT ( 1HO , ' 6. ANALYTICAL MODE-SHAPE VECTORS ', 

A ' (GRID PT / DI S PLACEMENT ) ' ) 

2000 FORMATdHO, 'ANALYTICAL MODE', 14, 

A 200(/5X,5(I9,'-‘ .11,'/' ,1PE10.3) ) ) 

3000 FORMATdHO, '7. EXPERIMENTAL MODE-SHAPE VECTORS ', 

A ' (GRID PT/DISPLACEMENT) ' ) 

4000 FORMATdHO, 'EXPERIMENTAL MODE', 14, 

A 200(/5X,5(I9, '-' ,I1,'/',1PE10.3>)> 

5000 FORMAT (1H1 , 'NASTRAN MODAL ANALYSIS - MODAL SURVEY STATISTICAL 
A 'CORRELATION', /1X,A9> 

END 
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SUBROUTINE RM8DIFF < AVECT , EVECT , NDOFI , RMS A ,RMSE, CORREL , 

A S,DIFF, DIFABS) 

C 

C CALCULATE RELATIVE (TO RMS) DIFFERENCE BETWEEN ANALYTICAL AND 
C EXPERIMENTAL DISPLACEMENT FOR EACH GRID POINT. 

C CALCULATE ABSOLUTE DIFFERENCE SCALED TO ANALYTICAL, 

C 

REAL AVECT(A) ,EVECT( A) ,DIFF(A) ,DIFABS(A) 


C 


IF (NDOFI .EQ. 0) RETURN 
DO 100 I»1 , NDOFI 
IF (CORREL .GE. 0.0) THEN 
DIFF(I) n 0.0 

IF ( (RMSA ,NE. 0.0) .AND, (RMSE ,NE. 0,0) ) 

A DIFF(I) = AVECT ( I ) /RMSA - EVECTf I) /RMSE 

DIFABS(I) =0.0 
IF (S .NE. 0.0) 

A DIFABS(I) a ( AVECT ( I ) - EVECT(I)) / S 

ELSE IF (CORREL .LT. 0,0) THEN 
C CORRECT FOR 1 BO-DEGREE PHASE SHIFT 

DIFFCI) « 0.0 

IF ( (RMSA .NE. 0.0) .AND. (RMSE .NE. 0.0) ) 

A DIFF(I) = AVECT { I ) /RMSA + EVECT( I ) /RMSE 


DIFABS ( I ) » 0.0 ;j 

IF (S .NE. 0.0) il 

A DIFABS ( I ) = (AVECT(I) + EVECT(I>) / S ; 

END IF ! 

100 CONTINUE 

RETURN .' 

END f 


j 
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SUBROUTINE SWAP(I,J) 

C 

C INTERCHANGE THE VALUES IN TWO 4-BYTE VARIABLES 

C 

ITEMP = I 
I = J 
J ■ ITEMP 
RETURN 
END 
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